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t  .0  Spacecraft  Charoin 


Magnetospher i c  substorms  at  high  altitudes  are  characterized  by 
hot  tenuous  plasmas  which  can  charge  a  satellite  negatively  to 
many  kilovolts.  The  reliable  estimation  of  the  charging  effects 
on  present  and  future  space  vehicles  has  been  a  concern  of  the 
Air  Force  and  NASA  for  several  years,  and  a  subject  of  consid¬ 
erable  RftcD  activity.  The  inclusion  of  all  factors  that  are 
significant  for  spacecraft  charging  presents  a  formidable 
problem.  Complex  three  dimensional  geometries  with  iiverse 
surface  materials  and  internal  electrical  pr open t i e  nav'e  to  be 
taken  into  account.  The  surface-plasma  interactio  jhich 
involve  secondary  emission,  backseat  ter i ng ,  photoer  sion  and 
charge  transport  must  be  considered  in  detail .  Th  gneto- 
static  environment  may  include  solar,  magnetic  field,  substorm 
and  space  charge  effects.  Charging  dynamics  can  exist  due  to 
varying  photoemission  over  a  rotating  or  eclipsing  satellite, 
changing  plasma  conditions,  active  control  by  electron  or  ion 
beams,  and  switching  of  electrically  connected  satellite 
configurations.  Investigators  at  AFGL  are  widely  involved  in 
SCATHA  satellite  and  other  plasma  physics  research  experiments. 

For  a  number  of  years  the  analytical  tool  which  most  compre¬ 
hensively  represents  these  processes  has  been  the  NASCAR  program 

(14) 

developed  by  S-CUBED.  ’  This  program  represents  an  attempt 

to  realistically  simulate  and  depict  spacecraft  charging  by 
numerically  solving  the  full  three  dimensional  problem.  The 
NASCAR  model  uses  a  quasi-static  Vlasov  equation  approximation  to 
obtain  surface  charging  fluxes,  given  potential  distributions  and 
the  ambient  environment  conditions.  Following  calculation  of 
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dielectric  charge  transport  and  surface  charge  distributions,  the 
spatial  potential  problem  (Laplace's  equation)  is  solved  using  a 
finite  element,  nested-mesh  technique.  The  iterative  cycle  then 
returns  to  the  Vlasov  flux  charging  equation.  Elaborate  plasma 
and  material  characteristics  may  be  formulated  in  NASCAP,  along 
with  complex  geometries  and  electrical  connectivities  of  the 
surfaces . 

Thus  NASCAP  is  a  major  evolving  engineering  level  program,  as 
well  as  a  tool  for  research  on  plasma  sheath  models  and 
spacecraft  charging  experiments  such  as  SCATHA.  The  NASCAP 
program  was  installed  at  AFGL  in  1978  for  the  purpose  of 
verification  and  validation,  and  as  a  means  of  support  for 
investigations  of  charging  phenomena. 

Since  installation,  the  program  has  been  run  extensively  to 
analyze  charging  effects  and  to  carry  out  simulation  studies  for 
SCATHA  and  other  space  vehicles.  The  early  worK  focused  on  the 
identification  of  program  bugs  and  computational  instabilities. 
These  studies  were  served  by  examining  the  charging  behavior  of 
simple  objects  such  as  NASCAP  quasi  spheres  and  cyl’nders. 

Further  testing  of  the  NASCAP  program  was  done  by  calculating 

sunlight  charging  of  quas 1 -spheres  for  fixed  sun  direction  and 

( 3) 

rotational  motion.  Preliminary  SCATHA  simulations  were 

^  0  '^ 

completed  assuming  simplified  plasma  environments.  *  These 

runs  included  quasi  spheres ,  2-grid  and  4-gnd  SCATHA  models. 

The  effect  of  electron  emitter  beams  on  SCATHA  charging  was  also 
examined.  '  The  most  recent  SCATHA  simulations  have 
implemented  time  dependent  specification  of  the  plasma 
environment  and  of  solar  flux.  ’ 
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1.1  Operational  System  and  Simulation  Procedures 

The  NASCAP  source  code  is  provided  by  S-CUBED  on  magnetic  tape  in 
blocked  EBCDIC  -format.  A  number  o-f  steps  are  necessary  in  order 
to  set  up  the  system  at  AFGL : 

a>  read  the  source  tape  and  convert  to  CDC  compatible  -format 

b)  convert  the  NASCAP  routines  to  UPDATE  format  and  compile 

c)  form  a  segmented  version  of  the  program  which  is  small  enough 
to  fit  on  the  CDC  6600 

The  final  segmented  binary  file  requires  276  KQ  to  run.  Also 
contained  on  the  source  tape  are  three  NASCAP  auxiliary  programs: 

PLOTREAD  -  reads  plot  files  created  by  NASCAP  and  outputs 
mi crof i che 

TERMTALK  -  reads  the  results  files  created  by  NASCAP,  giving 
access  to  summary  information. 

MATCHG  -  uses  NASCAP  algorithms  to  calculate  charging  of 
material  surfaces. 

In  addition,  S-CUBED  provides  an  auxiliary  file,  NASFILES,  which 
contains  default  values  for  the  surface  material  parameters,  and 
object  definition  specifications  for  SCAT HA  models. 

A  number  of  major  revisions  of  the  NASCAP  code  and  its  auxiliary 
programs  have  been  installed  on  the  CDC-6600  system. 

A  comprehensive  system  has  been  implemented  to  facilitate  program 
maintenance  and  operational  logistics.  Considerations  provided 
for  at  AFGL  are: 
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a>  coordination  and  compatibility  with  the  UNIUAC  1180  version 
of  NASCAP  at  S-CUBED) 

b)  prompt  revision  of  source,  relocatable,  and  absolute  codes 

c)  tape  storage  of  source  codes  and  results  files 

d)  generation  of  microfiche  plots  using  the  auxiliary  program 
PLOTREAD 

e)  maintenance  and  modification  of  auxiliary  programs  MATCHG 
and  TERMTALK. 

f)  availability  of  SUATEK  plotting  routines  for  display  of 
MAT CHG  and  TERMTALK  results 

g)  maintenance  and  modification  of  program  ALLES  which  calculates 
Maxwellian  fits  to  spectra,  for  input  to  NASCAP 

h)  generation  of  tapes  and  associated  materials  for  external 
distribution  of  programs 

The  NASCAP  simulation  of  a  complex  charging  event  entails 
considerable  operational  management.  Input  files  containing  the 
material  parameters  and  the  object  definition  information  are 
required.  A  flux  definition  file,  which  describes  the  plasma 
environment,  must  be  constructed.  Provision  for  inputting  the 
flux  files  to  the  auxiliary  MATCHG  program  for  trial  run  testing 
is  advantageous.  It  is  necessary  to  form  a  run  options  file 
from  the  NASCaP  Keyword  list,  in  accordance  with  the  rules  o+ 
syntax  and  precedence.  The  run  options  file  may  contain 
repetitive  modules  for  generating  variable  sunlight  intensity  or 
rotational  motion.  Special  card  blocks  may  be  needed  for 
specification  of  onboard  beam  emitter  or  detector  character¬ 
istics.  After  completion  of  a  NASCAP  run  the  printed  output  and 
microfiche  must  be  examined  and  interpreted,  and  the  restart 
files  either  archived  or  reset  for  a  continuation  run. 

Inspection  of  the  restart  file  information  and  plotting  of 
summary  and  comparative  results  is  made  available  through  the 
auxiliary  TERMTALK  program.  Figure  1.1  shows  a  schematic  of  the 
AFGL  operational  system  for  a  NASCAP  simulation  study. 
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Figure  1.1  NASCA?  Simulation  Operational  Syataa 

15 


In  the  SCATHA  simulations,  the  ambient  plasma  has  been 
represented  by  single  or  double  Maxwellian  -fits  to  measured  data, 
as  determined  by  the  ALLES  program.  The  material  properties 
have  been  set  to  de-fault  values  provided  by  S-CUBED.  For  the 
vehicle  object  specification,  the  1,  2,  3  and  4-grid  models  of 
the  SCATHA  satellite  have  all  been  utilized. 

The  1-grid  and  3-grid  models  of  SCATHA  are  shown  in  Fig.  1.2. 

The  main  difference  between  these  two  models  is  that  the  3-grid 
version  has  a  more  detailed  representation  of  the  external  booms. 
Since  the  computation  time  goes  up  repidly  with  grid  size,  most 
of  the  analysis  has  been  done  with  the  1-grid  model. 


1.2  Plasma  Environment  Specification 

The  accurate  simulation  of  charging  events  requires  a  detailed 
representation  of  the  ambient  plasma.  This  is,  in  itself,  a 
difficult  problem  and  a  substantial  effort  has  been  applied  to 
the  development  of  a  program,  ALLES,  to  generate  parameter  files 
to  specify  the  environment.  The  program  makes  double  Maxwellian 
fits  to  measured  spectra  data,  with  the  plasma  temperatures  and 
densities  as  the  floating  parameters.  . 

Plots  of  the  individual  fits  versus  the  data,  and  summary  plots 
of  the  fit  parameters  versus  time  are  available.  The  following 

describes  the  fitting  procedures  that  were  used  in  the  analysis 

<9) 

of  the  day  114  <1979)  charging  event, 
obtained  from  the  SC9  experiment. 


based  on  spectra 


The  SC9  plasma  data  resides  on  magnetic  tape  in  the  -form  of 

electron  and  ion  detector  counts.  The  data  was  recorded  at  16 

sec  intervals  and  the  counts  are  distributed  into  64  energy  bins 

ranging  -from  0  to  80  KeM.  In  ALLES  the  raw  counts  are  converted 

3  6 

to  phase  space  distribution  ■functions  f  (units  are  sec"/cm  )  by 
multiplying  by  appropriate  calibration  and  phase  space  factors. 
Figure  1.3  shows  a  typical  spectrum  <UT=27368  secs).  The  lower 
two  plots  give  ln(f)  for  the  electrons  (left)  and  ions  (right). 
The  upper  plots  give  expanded  views  of  the  low  energy  region  (0- 
10Ke<v»>  . 

In  the  analysis,  the  15  lowest  energy  bins  are  disregarded. 

Also,  the  data  points  are  corrected  for  the  shift  in  the  spectra 
due  to  vehicle  charging.  In  the  plots,  the  circles  show  the 
shifted  points  and  the  crosses  denote  the  unshifted  distribution. 
The  method  employed  to  fit  the  shifted  spectra  is  outlined  below: 

1)  The  spectra  are  first  split  int  high  and  low  regions  (the 
1 ow  energy  range  extends  from  vehicle  ground  plus  2  KeV 
up  to  the  separation  cut  which  was  set  at  15  KeO) . 

2)  The  high  energy  range  (  >15  KeV  )  is  fitted  with  a  single 
Maxwellian  function. 

3)  This  high  energy  fit  is  extrapolated  back  into  the  low 
region  and  subtracted  from  the  data. 

4)  A  second  Maxwellian  is  fitted  to  the  difference  spectrum 
in  the  low  energy  region. 
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A  virtue  of  the  above  method  is  that  it  clearly  separates  the 
fits  into  high  and  low  energy  components.  However,  when  the 
method  was  applied,  it  was  found  the  the  solutions  reached  could 
be  very  sensitive  to  the  location  of  the  separation  cut. 
Moreover,  sometimes  no  reasonable  physical  values  were  obtained 
from  the  fit.  To  remedy  this  situation,  a  second  level  of 
analysis  was  employed.  The  results  from  the  above  method  were 
input  as  starting  values  for  a  general  least  squares  fitting 
routine  (if  the  first  method  fails,  the  values  from  the  last 
previous  successful  case  are  input).  This  additional  least 
squares  fitting  was  invoked  whenever  it  could  improve  the  RMS 
error  of  the  fit. 


A  typical  fit  to  the  data  is  displayed  by  the  solid  lines  in  the 
lower  plots  of  Figure  1.3.  The  numbers  in  the  upper  corners 
give  the  fit  parameters  (temperatures  and  densities  corresponding 
to  the  low  and  high  regions).  The  solid  lines  in  the  upper 

plots  show  low  energy  single  Maxwellian  fits  for  comparison  which 
were  not  used  in  the  NASCAP  simulation. 


In  Figures  1.4  (a) -(d)  we  summarize  the  derived  plasma  parameters 
for  the  day  114  full  eclipse  period.  Although  there  are  some 
fluctuations,  the  temperatures  and  densities  are  seen  to  be 
relatively  flat  over  this  time  period. 


Before  being  used  for  NASCAP  simulations,  a  further 
made  to  the  double  Maxwellian  fit  parameters.  The 
densities  were  renormalized  to  the  ion  densities  to 
charge  neutralization  and  a  correction  was  made  for 
content  of  the  plasma. 


change  was 
el  ec  tr  on 
provide  for 
the  02 
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i .  3  Materials  Charging  Studies 


The  auxiliary  program  MATCHG  provides  an  effective  instrument  for 
studying  the  interactions  between  a  vehicle  surface  and  a  given 
plasma  environment.  The  code  calculates  the  equilibrium 
potential  for  any  material,  from  the  current  balance  equation, 
for  cases  where  surface  conductivity,  photoemission,  and  multi¬ 
dimensional  effects  may  be  neglected.  The  MATCHG  program  has 
been  run  extensively  to  investigate  the  eclipse  charging  behavior 
of  many  different  ma t er i al s . *  The  code  has  been  utilized  for 

sensitivity  studies  of  the  material  properties,  comparisons  of 
the  potentials  for  single  and  double  Maxwellians,  testing  of 
alternate  formulations  for  secondary  emission  equations,  and 
various  other  studies  of  the  surface-plasma  interaction.  The 
MATCHG  code  has  been  adapted  so  that  it  can  read  in  the  NASCAP 
flux  definition  files.  In  this  way,  the  ALLES  fits  to  spectra 
data  can  be  input  to  MATCHG  and  time  histories  of  the  eclipse 
charging  voltages  can  be  produced.  The  results  can  be  used  to 
compare  the  effects  of  different  fitting  methods  on  the  charging 
strength,  or  for  making  trial  passes  through  the  data,  before 
employing  the  NASCAP  code  which  consumes  much  more  computer  time. 
Figure  1.5  shows  a  time  history  of  the  MATCHG  equilibrium 
voltages,  corresponding  to  ALLES  tits  to  the  day  113,  1981 
spectra  data.^1  *  Time  histories  can  also  be  generated  for  all 
of  the  contributing  currents  at  the  material  surface  (primary 
plasma  currents,  secondary  emission  and  backscatter  currents,  and 
bulk  conductivity  currents).  This  information  can  be  quickly 
obtained  and  gives  an  indication  of  the  relative  importance  of 
these  effects  throughout  the  charging  simulation. 

A  series  of  MATCHG  runs  have  been  done  to  investigate  systemati¬ 
cally  the  threshold  conditions  for  materials  charging  based  on 

<  14) 

the  current  balance  condition.  The  dynamics  of  the  charging 
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Figure  1,5  MATCHG  Potentials  versus  UT  for 

Single  Maxwellian  Fits  (Day  113/81) 


is  illustrated  in  Figure  1.6(a)  which  shows  the  MATCH6  electron 

and  ion  currents  -J  (V)  and  J  (V)  versus  V  for  ASTROQ  material . 

e  i 

Each  curve  corresponds  to  a  different  Maxwellian  temperature 

<T  ,  T.  =  5,10,  15,  20,  25).  The  plasma  densities  are  taken  as 

n  “n^l.  The  intersection  of  any  two  curves  determines  the 

equilibrium  vol taqe  U.  The  curve  for  -J  (O)  with  T  =5  lies 

e  e 

below  the  origin  and  hence  can  never  intercept  a  J^V)  line. 

This  is  an  example  of  a  case  where  the  electron  temperature  is 
below  the  threshold  and  no  charging  can  occur.  The  remaining 
-Je<U>  curves  start  above  the  J.(U)  lines  at  0=0  and  charging 
takes  place. 

We  have  examined  the  threshold  conditions  for  charging  of  a 
surface  material  immersed  in  a  Maxwellian  plasma.  The  threshold 
relations  depend  on  the  plasma  temperatures  and  densities  and  on 
the  material  properties  through  the  secondary  emission  coeffi¬ 
cients  &^(T  ),  £®(T.)  and  the  backscatter  ratio  £B(T  >. 

e  i  e 

For  a  single  Maxwellian  plasma,  the  necessary  condition  for 
charging  is  that  the  electron  temperature  must  be  greater  than 
the  threshold  temperature  defined  by 

&<T  )  =  -1  +  fcS(T  )  «•  68<T  )  =  0 

e  e  e 

The  values  of  this  threshold  function,  calculated  by  MATCHG,  are 

shown  in  Figure  1.6(b)  for  five  different  NASCAP  materials. 

For  a  mu  1 t i -Maxwel 1 l an  plasma  at  least  one  of  the  electron 

temperatures  must  be  above  threshold  and  the  net  electron  induced 

current  must  be  negative.  To  arrive  at  necessary  and  sufficient 

conditions  for  charging,  one  has  to  take  into  account  the  ions  as 

well  as  the  electrons.  Explicit  formulas  have  been  given  for 

single  and  double  Maxwel lians.  For  a  double  Maxwellian,  the 

threshold  condition  defines  a  linearly  bounded  region  in  the 

density  space  <n  n  >  . 

el  e2 
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AND  ION  CURRENTS 


1.4  Special  Simulations 


The  NASCAP  model  has  been  applied  to  several  different  space 
vehicles.  A  simple  geometry  was  con-figured  in  order  to  repre¬ 
sent  the  Helios  1  satellite  and  its  external  antenna  boom.^®'> 

The  NASCAP  configuration  that  was  used  is  shown  in  Figure  1.7. 

The  satellite  body  is  a  2x2  mesh  cube.  Attached  is  a  single  16- 
segmented  cylindrical  boom  (radius  =  7mm)  which  is  connected  to  a 
separate  conductor  (other  than  the  vehicle  ground)  through  a  286 
pf  condenser  (the  Helios  antenna  capacitance).  In  order  to 
scale  the  system  to  correspond  to  the  Helios  dimensions,  the  mesh 
size  was  set  to  1  meter. 

Figure  1.7  also  shows  the  boom  potentials  calculated  during  one 
cycle  of  the  NASCAP  rotate  mode.  Two  curves  are  given:  one  for  a 
cell  at  the  base  of  the  boom,  and  one  for  a  cell  at  the  boom  tip. 
As  can  be  seen,  there  is  only  a  small  variation  of  potential 
volt)  along  the  boom.  Both  curves  show  the  expected  spin 
modulation.  When  the  boom  is  perpendi cul ar  to  the  sun 
direction  (♦  =  +90,-98°),  the  photo  current  is  a  maximum  and  the 
voltage  goes  positive.  When  the  boom  is  aligned  with  the  sun 
<♦  =  0,-188°)  the  cells  go  a  few  volts  negative. 

A  comparison  of  Figure  1.7  with  the  Helios  measurements  of  boom 
potential  indicate  that  the  amplitude  of  the  oscillations  are  of 
the  right  order  of  magnitude.  Considering  the  very  simplified 
model  that  has  been  employed  here,  the  above  agreement  is 
satisfactory.  It  demonstrates  that  the  NASCAP  program  can 
successfully  predict  the  spin  modulated  potential  variations  on 
a  rotating  boom,  and  that  the  details  of  the  model  are  not 
significant  for  this  case. 
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Figure  1.7 


Helios  Model  and  Boom  Potentials  vs  Rotation  Angle 


The  NASCAP  program  has  been  used  to  configure  a  model  of  the  GEOS 

<  13) 

spacecraft  and  SPA  sensor  unit.  The  GEOS  satellite  consists 

of  a  cylindrical  body  with  projecting  booms.  The  Suprathermal 
Plasma  Analyzer  detector  is  mounted  on  one  of  the  booms. 

The  corresponding  NASCAP  model  is  described  in  Figure  1.8(a). 

The  main  features  of  the  model  are: 

a)  The  satellite  body  is  built  up  from  5  vertically  stacked 
8x4  mesh  octagons.  The  mesh  size  was  taken  as  .2  m  so  the 
effective  diameter  is  8x.2  =  1.6  m  and  the  height  is 

5x . 2  =  1.0m. 

b)  Two  external  cylindrical  booms  of  length  llx.2  =  2.2  m 
extend  out  from  the  top  of  the  satellite  body. 

c)  The  SPA  sensor  unit  is  represented  by  a  Ixlxl  mesh  cube. 

The  effective  dimensions  are  thus  20x20x28  cm,  which  is 
slightly  larger  than  actual  size.  The  SPA  cube  was  put  on 
a  separate  conductor  to  allow  for  independent  biasing. 

d>  The  1 ower  3  mesh  units  <.6m>  of  the  satellite  body  is 
covered  with  the  material  Index. 

e)  All  other  surfaces  (the  top  of  the  satellite,  the  booms,  the 
SPA  unit)  are  covered  with  Goldpd. 

Sunlight  charging  of  a  GEOS  satellite  model  for  three  repres¬ 
entative  voltages  on  the  SPA  detector  were  studied. 


28 


!*■>/  •  I  HA<  •  M  «  V  i 

ll'M  H|M  >  ‘  ‘I  "'A  •  -l« 


ft  -  ' i  <N  •  ■  S'  A  *  F 


Figure  1.8(b)  GEOS  Potential  Contours  (Top  View) 
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The  primary  questions  to  be  examined  were: 


1)  Are  the  local  electrostatic  -fields  around  the  detector 
significantly  perturbed  by  changes  in  sensor  bias? 

2)  Are  the  trajectories  o-f  low  energy  particles  near  the  SPA 
instrument  significantly  altered  by  the  sensor  bias? 

When  the  model  is  exposed  to  the  expected  plasma  environment  and 
sunlight,  the  vehicle  frame  potential  goes  a  few  volts  positive. 
If  the  bias  on  the  SPA  surface  is  varied  from  0  to  -28  volts,  the 
local  field  near  the  unit  is  modified  and  there  is  a  significant 
change  in  the  low  energy  trajectories.  Figure  1.8(b)  shows  a 
top  view  of  the  potential  contours  around  the  SEOS  model.  The 
influence  of  the  detector  biasing  on  the  contours  is  clearly 
evi den  t . 


The  proposed  Space  Based  Radar  structure  was  modeled  by 
(12) 

NASCAP.  The  SBR  vehicle  is  composed  of  three  main 

components: 

1)  the  lower  side  package  (LSP)  consisting  mainly  of  the  large 
antenna  array. 

2)  the  upper  side  package  (USP)  consisting  of  two  gimbaled  solar 
arrays . 

3)  the  mast  connecting  LSP  to  USP. 
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The  corresponding  NASCAP  model  is  given  in  Figure  1.9(a). 

The  mesh  size  is  8  meters.  The  large  antenna  <LSP>  is 
represented  by  an  octagon  of  diameter  8  meshes  =  64  meters 
(actual  diameter  =  71  meters). 

The  mast  is  a  boom  cylinder  of  diameter  2  meters  and  length  12 
meshes  =  96  meters  (actual  length  182.5  meters).  The  solar 
arrays  ( USP)  are  represented  by  a  2x2  mesh  =  16  meter  x  16  meter 
sauar e . 

The  surface  material  on  the  large  antenna  (LSP)  is  mainly  KAPTON 
with  a  small  amount  of  ALUMINUM  edge  tape  and  the  USP  solar  array 
is  assumed  to  be  borosilicate  g1ass.(1>  For  the  NASCAP  model  we 
have  assumed  that  the  LSP  and  the  mast  are  covered  with  KAPTON 
and  the  USP  with  SOLAR  material . 

The  SBR  will  be  in  polar  orbit  going  through  the  auroral  region 
at  altitude  %10,000  Km.  The  LSP  antenna  always  points  earth¬ 
ward,  while  the  gimbaled  solar  arrays  align  themselves  with  the 
sun  direction.  Relative  to  a  coordinate  system  fixed  on  the  SBR 
the  sun  appears  to  rotate  through  360°  on  each  orbit.  We  chose 
two  representative  positions,  at  +45°  and  -45°,  for  the  NASCAP 
simulations. 

Two  test  plasma  cases  were  set  up: 

1)  a  "cool"  plasma  with  T  =5,  T.=18 

e  l 

2)  a  "hot"  plasma  with  Te=18,  T^=28 

-3 

The  densities  were  assumed  to  be  1  cm  ,  In  the  cold  plasma  the 
materials  ALUMINUM  and  KAPTON  would  charge  (in  eclipse)  but  SOLAR 
would  not.  For  the  hot  plasma  all  the  materials  would  charge. 
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The  NASCAP  runs  that  were  considered  are  summarized  in  Table  1.1. 
The  vehicle  ground  voltage  and  the  potentials  on  selected  surface 
cel  1 s  are  listed. 

For  case  ( 1) ,  the  cool  plasma  is  used  and  the  sun  is  at  angle 
♦45°.  The  vehicle  ground  charges  to  -1.8  KU .  The  potential 
contour  plots  indicate  that  charge  builds  up  on  the  shade  side  of 
the  LSP .  The  solar  array  does  not  charge  and  therefore  the  field 
lines  are  unperturbed  near  the  USP. 

Case  <2)  is  the  same  as  case  < 1)  except  that  the  hot  plasma 
replaces  the  cool  plasma.  The  vehicle  ground  thus  increases  to 
-5.1  KU .  The  potential  contours  are  similar  to  those  of  case  (  1) 
but  now  the  USP  charges  up  and  perturbs  the  potential  lines  in 
its  vicinity. 

Case  <3)  differs  from  case  <2)  in  that  the  sun  angle  is  changed 
from  +45°  to  -45°.  Here  the  contour  lines  build  up  on  the 
topside  of  the  LSP,  which  is  now  in  the  shade.  The  vehicle 
ground  is  about  the  same  as  for  case  (2) . 

The  simplified  NASCAP  model  of  the  SBR  provides  a  possible 
scenario  for  vehicle  charging.  As  it  follows  a  polar  orbit, 
the  top  and  bottom  sides  of  the  LSP  revolve  into  and  out  of  the 
sunlight.  If  there  is  a  low  density  ambient  plasma  present, 
with  electron  temperature  greater  than  3KeU ,  then  the  shade  side 
of  the  LSP  can  charge  up  to  the  equivalent  eclipse  equilibrium 
voltage.  If  the  plasma  temperature  is  high  enough  (%6KeV>,  then 
the  shade  side  of  the  USP  will  charge.  The  potential  contours 
are  of  the  type  shown  in  Figure  1.9(b). 
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TABLE  1.1  NASCAP  RESULTS  FOR  SBR  CHARGING 


CASE  (  l>  I  CASE  (2)  )  CASE  (3) 


PLASMA 
SUN  ANGLE 
ITERATION  STEP 
TIME 


1 . 8  2x  10' 


3 . 2x  10 


3.2x  10 


GROUND  VOLTAGE  -1 


LARGE  ANTENNA 
LSP,  bottom 


LARGE  ANTENNA 
LSP,  top 


MAST-  bottom 
mi  ddl e 
top 


SOLAR  ARRAY 
USP,  bottom 


SOLAR  ARRAY 
USP,  top 


C-ELL# 

MATERIAL 

V 

V 

-1 

-1S03 

-514? 

2 

KAPTON 

-2681 

-7467 

1 

H 

-1329 

-3802 

161 

II 

-1006 

-37?5 

166 

u 

-735 

-22?8 

172 

u 

-473 

-2638 

65 

SOLAR 

-48? 

-3087 

64 

SOLAR 

-408 

- 1?22 

3 


SCATHA  charging  in  sunlight  has  been  studied  -for  a  "worst  case" 
plasma  environment,  based  on  SC5  measurements.^'  1  In  the 

NASCAP  simulations  an  averaged  double  Maxwellian  -fit  to  the 
plasma  spectrum  was  used.  Due  to  uncertainties  in  the  envi¬ 
ronment  parameters  we  have  not  attempted  to  predict  absolute 
potentials  with  NASCAP  but,  instead,  have  looked  at  the  behavior 
of  SCATHA  charging,  under  a  variety  of  input  conditions,  tor 
comparative  studies.  The  results  ot  the  NASCAP  runs  are  given 
in  Tables  1.2  and  1.3.  The  cases  that  have  been  considered  are 
labelled  from  1  to  8  as  indicated. 

Cases  < 1)  and  (2)  compare  eclipse  charging  tor  the  1-grid  SCATHA 
model  using  secondaries  ot  normal  incidence  (NORMAL)  or  with  an 
isotropic  distribution  WANGLE).  It  can  be  seen  that  the  ground 
voltage  is  decreased  trom  -10.9  KU  to  -7.0  KV  when  ANGLE  replaces 
NORMAL.  This  reduction  is  due  to  the  tact  that  more  secondaries 
are  emitted  when  the  particles  impinge  on  a  surface  isotropi¬ 
cally.  The  ANGLE  result  is  in  good  agreement  with  the  voltage 
that  was  actually  measured  -8  KLD  in  the  eclipse  that  occurred 
about  25  minutes  later. 

Cases  (3)  and  (4)  again  make  the  comparison  between  NORMAL  and 
ANGLE  secondaries,  but  for  a  sunlight  environment.  Here  the  SPIN 
mode  ot  NASCAP,  set  to  give  spin  averaged  illumination  around  the 
X-axis,  has  been  employed.  Again  there  is  a  large  difference 
(-3.2  to  -.74  KU)  between  NORMAL  and  ANGLE  results  for  the  ground 
voltage.  The  voltage  for  the  ANGLE  calculation  is  about  a  factor 
of  two  greater  than  than  the  maximum  measured  sunlight  potential 
< - . 34  KV>  . 

In  the  preceding  runs,  it  was  assumed  that  there  was  no  biasing 
of  the  conductors  2  and  3  of  the  1-grid  model.  To  be  consistent 
with  the  S-CUBED  configuration,  these  conductors  were  biased  to 
zero  volts  relative  to  conductor  1.  Case  (5)  mss  run  to  test 
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TABLE  1.2  SUMMARY  CASES  i  to  4  (WORST  CASE.) 

ECLIPSE  SPIN 


r 

T -  T 

1 

[ 

- r 

j 

\ 

L 

CELL 

IMATERIAL 

1 

i 

NORMAL  1 

1 

i 

ANGLE  I 

1 

i 

NORMAL 

j 

ANGLE  | 

i 

r 

-1 

! 

i  GROUND 

1 

1 

-10 ,886 l 

! 

i 

-7039  I 

1 

-3155 

i 

i 

-736  | 

I 

-2 

1 

1  BOCiMT  IPS 

1 

1 

-14,473  I 

I 

-9885  1 
! 

-1496 

-498  | 

1 

-4 

1 

IREF  BAND 

1 

1 

-11,512! 

1 

1 

-7559  I 

1 

-2576 

1 

-1114  | 

! 

279 

1 

1  KAPTON 

1 

-14,450  I 

-9891  I 

1 

-2122 

-614  | 

1 

SC  1  -  1  1 

31  1 

1 

1  SOLAR 

1 

! 

-15, 160  I 

I 

-9731  ! 

1 

-2136 

1 

-6,2  | 

309 

1  INDOX 

1 

-10 ,890  i 

! 

-7039  ! 
i 

-3155 

-736  | 

419 

1 

I  GOLD 

1 

! 

-1  1 ,860  I 

1 

-7695  I 

1 

-2108 

1 

-612  I 

1 

SC  1-2  1 

135 

I  KAPTON 

1 

1 

-14,450  1. 

-9891  *  I 

1 

-2174 

-802  I 
1 

352 

1 

1  GOLD 

I 

-10 ,890  | 

1 

-7039  I 

1 

-3155 

1 

-736  I 

1 

SC  1-3  i 

353 

I 

1  TEFLON 

1 

1 

-13,250  I 

| 

1 

-8168  1 

1 

j 

-13,220 

-8145  | 

| 

362 

|  KAPTON 

1 

-14,450  1 

-9891  I 

1 

-14,440 

-9865  I 

1 

363 

1  ASTROQ 

1 

-10  ,930 1 

1 

-7060  1 

1 

-3455 

1 

-1081  I 

1 

BOTTOM  I 

2 

1 

1  WHITEN 

1 

1 

-11,5501 

1 

-7722  1 

1 

-5285 

I 

-2791  I 
1 

TOP  1 

359 

1 

1  GOLDPD 

1 

1 

1 

-10,890 1 

! 

-7039  I 

1 

1 

-3155 

I 

-736  I 
1 
i 

JOB  1 

1 

1 

1 

1 

1  1 

1  1 

TZWCEGR  ITZWCE2F  1 

1  1 

TZWCS13 
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TA8LE  1,3  SUMMARY  CASES  5  to  8  (  WORST  CASE) 
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the  effect  of  this  change.  By  comparing  cases  (4)  and  <5)  we 
see  that  this  modification  makes  negligible  difference  to  the 
charging.  On  subsequent  runs,  the  zero  biasing  of  conductors  2 
and  3  was  maintained. 

Figure  1.2  shows  the  NASCAP  1-grid  model  of  SCATHA  and  the 
corresponding  3-grid  model.  The  only  difference  between  the 
models  is  that  the  external  booms  have  been  modified.  In  the  3- 
grid  model  the  booms  extend  out  to  actual  length,  and  the  surface 
material  is  changed  from  BOOMAT  to  KAPTON.  The  cell  numbering 
differs  slightly  for  the  two  models.  In  column  1  of  Table  1.2 
both  numbers  are  given  < 1-grid  first),  where  differences  occur. 

The  3-grid  model  was  run  in  eclipse,  and  it  was  found  to  give 
virtually  the  same  potentials  as  for  the  1-grid  model  (case  2>  . 
This  result  is  expected  since  the  plasma  is  strong  enough  to 
charge  the  surface  cells  independently.  The  situation  can  be 
quite  different  in  a  weak  plasma  where  the  "boot  strap"  mecha¬ 
nism  (charging  cells  influencing  neighboring  non-charging  cells) 
can  lead  to  substantial  deviations  in  the  potentials  obtained 
from  the  two  models. 

Cases  (5)  and  ( 6)  compare  the  results  for  sunlight  charging  (in 
the  SPIN  mode)  with  the  l-gnd  and  3-grid  models.  It  is  evident 
that  the  3-grid  mode)  charges  higher  (-1.37  KV  versus  -.69  KV>  . 
This  result  holds  true  also  for  cases  (7)  and  (8)  which  make  the 
same  comparison  relative  to  the  ROTATE  mode  (-1.9  KV  versus  -1.23 
KV) .  The  potential  contours  for  the  3-grid  model  are  quite 
different  in  the  region  of  the  booms  from  those  of  the  l-gnd 
model.  Otherwise,  the  main  features  of  the  charging  are  the 
same  i .e.  the  WHITEN  cells  in  the  bottom  shade  and  the  SC1-3 
patch  in  the  top-side  shade  are  at  the  highest  potentials. 
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By  comparing  cases  (5)  and  < 7)  to  cases  (6>  and  ( 8>  we  see  that 
the  vehicle  ground  tends  to  be  higher  -for  the  rotate  mode  than 
with  SPIN  -for  both  the  1-grid  mode)  (-1.23  KV  versus  -.69  KV.)  and 
the  3-grid  model  (-1.91  KV  versus  -1.37  KV) .  This  behavior  is 
consistent  with  earlier  NASCAP  studies'' which  gave  the  same 
result  -for  a  uniform  quasi  sphere. 

The  NASCAP  comparative  runs  illustrate  the  following  charact¬ 
eristics  of  SCAT HA  charging: 

1)  The  ANGLE  mode  for  secondary  emission  produces  less  charging 
than  NORMAL. 

2)  The  1-grid  model  charges  up  less  than  the  3-grid  model. 

3)  The  biasing  of  conductors  2  and  3  on  the  1-grid  model  has 
little  effect 

4)  The  SPIN  mode  tends  to  give  lower  charging  than  for  ROTATE. 

In  the  ROTATE  runs,  the  amplitude  of  the  spin  modulation  can 
depend  on  the  DVLIM  parameter  which  may  have  to  be  adjusted  from 
the  default  value. 

The  NASCAP  calculation  of  the  vehicle  ground  potential  in 
sunlight  is  ^-1306  V  (average  over  cases  (5)  to  <S)>,  which  is 
higher  than  the  maximum  measured  value  of  -340  V.  The  eclipse 
charging  result  of  »;  7.0  KV  (from  case(2>>  is  less  than  the 
estimated  value  of  -16  KV  (from  Table  1).  On  the  other  hand  the 
NASCAP  result  is  quite  close  to  the  actual  measured  voltage  of 
%  8  KV  for  the  eclipse  which  occurred  later.  The  SSPM  results 
for  sunlight  charging  are  mixed:  the  calculated  voltages  are 
found  to  be  high  for  60. .0  and  SOLAR,  low  for  QUART 2 ,  and  in  rough 
agreement  wi th  the  KAPTDN  and  TEFLON  samples. 
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1.5  Extended  Simulation  Studies 

The  P78-2  SCATHA  satellite  was  launched  in  January  1978,  into  a 
near-geosynchronous  orbit.  The  satellite  carries  thirteen 
separate  experimental  payloads  -for  analyzing  spacecraft  charging 
effects.  The  data  from  the  UCSD  charged  Particle  Experiment 
(SC9)  was  used  to  determine  periods  of  strong  eclipse  charging, 
by  looking  for  sharp  peaks  in  the  high  resolution  electrostatic 
analyser  (ESA)  ion  count  measurements.  The  channel  number  of 
the  peak  gives  an  estimation  of  the  shift  in  energy  due  to 
particle  acceleration  into  the  detector,  and  thus  of  the 
spacecraft  potential.  This  effect  is  clearl y  evident  in  the 
upper -right  plot  of  Figure  1.3  where  there  is  a  clear  peak  in 
the  ion  distribution  function  at  about  6.2  KeV  corresponding  to 
vehicle  charging  of  6.2  Kv.  Figure  1.10  gives  a  summary  of  the 
spacecraft  potential  versus  UT  for  the  most  active  charging  days 
of  the  1979  eclipse  seasons.  The  "worst  case"  charging  event 

occurred  on  day  114,  1979  and  this  was  chosen  for  a  NASCAP 

.  .  .  .  .  (9) 

simulation  study. 

On  day  114,  the  SCATHA  satellite  was  in  full  eclipse  from 
UT=26055  to  UT=28868.  The  NASCAP  simulation  runs  were  started 
at  UT=25950,  with  the  vehicle  initially  uncharged  and  the  sun 
intensity  equal  to  zero.  This  procedure  approximates  the  actual 
situation  where  the  vehicle  was  charged  slightly  (*•  466  volts)  in 
sunlight  and  the  sun  intensity  rapidly  decreased  to  zero. 

The  NASCAP  1-gnd  model  (see  Figure  1.2)  was  used  to  represent 
the  SCATHA  satellite.  The  charging  simulation  was  carried  out 
in  the  UPDATE  mode.  That  is,  on  each  time-step,  the  program 
searches  an  input  file  for  the  most  recent  spectrum  and  updates 
the  flux  description  to  the  corresponding  double  Maxwellian  fit. 
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Figure  1.11  displays  the  NASCAP  results  tor  the  complete  pass 
through  the  data.  As  can  be  seen  -from  the  plot,  the  overall 
agreement  with  the  measured  SC9  voltages  is  quite  good. 

The  potential  contours,  -for  a  point  mid-way  through  the  second 
pass  run,  are  shown  in  Figure  1.12.  Near  the  surface,  the 
contour  lines  tend  to  follow  the  satellite  geometry,  with  some 
deviations  due  to  the  booms.  Beyond  a  few  meters  from  the 
satellite,  the  monopole  field  dominates. 


The  most  complex  NASCAP  simulation  that  has  been  attempted  is  the 
day  113  (1981)  case,  where  transit  into  and  out  of  eclipse  was 
considered  /  Here,  in  addition  to  the  temporal  plasma 

environment,  we  have  time  dependent  sun  illumination  in  the 
pen  umbr  a  regions. 

The  SC9  instrument  sweeps  through  64  energy  channels  <0  <  E  <  82 
KeV) .  The  measurements  were  made  once  every  17  seconds  during 
the  day  113  eclipse  period.  The  ALLES  program  reads  the  SC9 
data  from  magnetic  tape,  converts  the  electron  and  ion  counts  to 
the  corresponding  phase  space  distribution  functions,  and  makes 
Maxwellian  fits  to  each  spectrum.  Figure  1.13  shows  the  SC9 
distribution  functions  versus  UT . 

To  obtain  an  overview  of  the  charging  strength  of  the  day  113 
plasma  environment,  a  number  of  trial  runs  were  done  with  the 
auxiliary  program  MATCHG.  The  MATCHG  program  was  adapted  so 
that  it  could  read  in  the  spectrum  fits,  one-by-one,  and 
calculate  the  equilibrium  potential  for  each  of  them.  The 
surface  material  was  chosen  to  be  SOLAR,  which  is  the  most 
prevalent  cell  on  the  NASCAP  model  of  SCATHA.  From  previous 
work  it  has  been  found  that  SOLAR  is  a  good  indicator  of  the 
eclipse  charging  behavior  of  the  1-grid  model  i.e.  when  SOLAR 
charges,  the  NASCAP  model  does  so  also,  usually  reaching  a  value 


30’/.  below  that  of  MATCHG. 


Day  114.  1979 
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Figure  1.13  Electron  and  Ion  Phase  Space 

Distribution  Functions  for  Day  113 
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These  studies  indicated  that  adequate  charging  voltages  in  the 
eclipse  region  could  be  obtained  by  using  the  single  Maxwellian 
tits  from  ALLES.  Figure  1.5  gives  the  MATCHG  equilibrium 
potentials  versus  UT  corresponding  to  these  tits.  In  order  to 
reduce  fluctuations,  the  electron  densities  and  temperatures  were 
smoothed  using  a  ten  point  running  average.  The  tinal  input 
temperatures  and  densities  are  given  in  Figure  1.14. 

To  conserve  computing  time,  the  NASCAR  1-grid  model  ot  SCATHA  was 
employed.  This  model  gives  an  adequate  representation  ot  the 
vehicle  geometry  and  surtace  materials  (the  material  properties 
were  assumed  to  have  the  detaul t  values  set  by  S-CUBED) . 

However,  the  booms  are  truncated  and  this  will  lead  to  some 
inaccuracy  in  the  sunlight  regions.  Also,  to  save  time,  the  sun 
illumination  was  represented  by  the  SPIN  approximation,  rather 
than  ROTATE.  Again  this  will  introduce  uncertainties  into  the 
calculation  of  sunlight  effects. 

The  NASCAR  simulations  were  carried  out  in  the  UPDATE  mode  i.e. 
on  each  iteration  the  program  searches  an  input  -file  containing 
the  spectra  -fits  and  picks  out  the  one  that  is  closest  (earlier.) 
to  the  current  time.  Also,  by  means  of  RDOPT  cards,  the  sun 
intensity  is  adjusted  on  each  step  to  define  the  variation  in  sun 
illumination  in  the  entrance  and  exit  penumbra  regions. 

The  NASCAR  time-step  was  set  to  17  secs,  which  is  the  interval 
between  SC9  measurements,  so  that  each  spectra  would  be  sampled 
at  least  once  during  the  run.  A  complete  pass  through  the  data 
--  into  and  out  of  eclipse  --  takes  about  6  hours  running  time  on 
the  CDC  6600,  using  the  1-grid  model . 

The  vehicle  frame  potential,  as  determined  from  SC9  measurements, 
is  plotted  versus  UT  in  Figure  1.15.  The  circles  represent  the 
measured  voltages  in  KM.  The  solid  line  in  the  figure  gives  the 
sunlight  intensity  during  passage  through  eclipse.  The  voltage 
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Figure  1.14  Corrected  and  Smoothed  Electron  Density 
and  Temperature  (Single  Maxwellian) 


plot  has  two  peaks;  the  first  peak  is  correlated  with  entry  to 
eclipse  and  the  second  one  occurs  in  the  exit  penumbra. 

A  straightforward  physical  picture  may  be  used  to  explain  the 
main  features  of  the  day  113  charging  event.  Assume  that  there 
was  present  a  severe  charging  plasma  environment.  In  sunlight, 
any  shaded  surfaces  would  charge  up,  but  those  cells  that  are 
exposed  to  the  solar  IA>  would  lose  photo-electrons,  so  that  the 
charging  would  be  inhibited.  The  net  effect  is  that  the  vehicle 
ground  typically  goes  a  few  hundred  volts  negative.  When  the 
spacecraft  enters  penumbra,  the  photoemission  is  reduced  and 
eventually  a  point  is  reached  where  the  plasma  current  dominates. 
The  vehicle  then  rapidly  charges  and  maintains  a  high  voltage 
through  eclipse.  In  the  exit  penumbra  region,  the  photoemission 
again  builds  up  and  the  charging  is  suppressed. 

Figure  1.15  gives  the  NASCAP  predictions  for  the  day  113  event 
assuming  a  DVLIM  parameter  of  580  <a  comparison  run  with  value 
5000  was  also  done).  The  following  aspects  of  the  simulations 
are  worthy  of  note: 

a)  in  sunlight,  the  vehicle  frame  potential  goes  a  few  hundred 
volts  negative,  as  expected. 

b)  when  the  sun  illumination  in  the  entrance  penumbra  is  suffic¬ 
iently  reduced  <  to  %  25'/.)  there  is  rapid  charging.  With  the 
smaller  CWLIM  setting  <598),  the  charging  peak  is  *  18KV,  in 
agreement  with  measurements. 

c)  in  the  full  eclipse  region,  the  vehicle  frame  falls  to  s?  6KV, 
which  is  a  few  kilovolts  higher  than  the  measurements.  This 
result  is  independent  of  DULIM. 
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d>  the  second  peak  structure  in  the  data  is  too  high  when 
DMLIM  =  5008,  and  slightly  low  -for  DVLIM  =  508. 

The  latter  case  should  be  more  accurate. 

e)  in  the  exit  penumbra,  the  solar  UV  brings  down  the  -frame 
potential,  roughly  in  accordance  with  measurements.  The 
region  is  almost  independent  of  IX>LIM  even  though  there  are 
frequent  cut-backs  with  DVLIM  =560  and  none  with  DO LIM  =5000. 


The  primary  weakness  in  the  simulation  appears  to  be  near  the 
entrance  penumbra  region  where  the  NASCAP  frame  potential  starts 
to  rise  later  than  for  the  measurements  and  then  goes  up  much 
more  rapidly.  This  is  an  inherently  difficult  part  of  the 
event  to  model,  requiring  precise  values  of  the  sun  illumination 
and  correlated  plasma  environment. 

The  SPIN  approximation,  which  averages  the  sun  intensity  over 
many  cycles  is  probably  inadequate  here  (the  jogs  in  the  measured 
voltage  during  charging  are  probably  due  to  rotational  effects). 
Also,  neglecting  the  low  energy  electrons  is  less  valid  in  this 
region  where  the  vehicle  potential  is  low.  For  these  reasons, 
it  is  not  surprising  that  the  NASCAP  code  does  not  track  the 
measurements  better  through  the  initial  charge-up  phase. 


The  overal 1  agreement  of  the  1-grid  frame  potential  with  the  SC9 
voltages  is  quite  good  considering  that: 

1)  we  have  neglected  the  low  energy  Maxwellian  components 

2)  there  are  uncertainties  in  the  materials  properties 
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3)  there  are  approximations  in  the  model  (truncated  booms, 

SPIN  mode  instead  of  ROTATE) 

4)  there  is  uncertainty  in  the  0-,  content  o-f  the  environment 

5)  there  were  large  -fluctuations  in  the  charging  strength  -from 
one  plasma  spectrum  to  the  next. 

The  first  point  above  is  the  most  crucial  one.  The  low  energy 
electrons,  in  particular,  form  a  critical  part  of  the  plasma 
specification.  But  the  low  energy  part  of  the  SC9  spectra  is 
not  Known  during  periods  of  strong  charging  because  these 
particles  are  repelled  by  the  electric  field  surrounding  the 
vehicle.  Thus  some  guess  has  to  be  made  as  to  the  shape  of  the 
low  energy  electron  spectrum.  If  we  use  extrapolated  fits  to 
the  region  (the  double  Maxwellian  fit)  then  we  find  that  the 
charging  is  suppressed.  If  we  neglect  the  1 ow  energy  components 
entirely  (the  single  Maxwellian)  then  we  get  charging  that  is  a 
few  kilovolts  too  high  in  the  central  eclipse  region.  A 
definitive  test  of  the  NASCAP  program  can  only  be  made  if  a 
complete  specification  of  the  particle  spectrum  is  available. 
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2  .  £'  I  on  i  z  a  t  i  on _ and  Charging  Effect-;.  of  Beam  Curr  en  t  s- 

When  a  charged  probe  is  placed  in  a  plasma  err  ■  1  f  onmen  *  .  i  t.  driH: 
■from  the  plasma  a  current  ct  oppos i  tely  charoed  particles. 

ince  particles  of  1  ike  charge  are  repel  led  bv  tne  pr ode, 
sheath  develops,  depleted  to  some  degree  i  n  the  repelled  species. 
I  he  application  ot  this  system  to  the  charging  of  beam  emi  Mine 

(  i  i 

eoacecr  at  t  was  -first  detailed  by  Beard  and  .Johnson  .  Later, 

the  physics  of  the  probe  sheath  was  studied  in  depth  bv  Lam. 

(  4 . 

More  recently,  Leadon ,  et.al .  developed  a  numerical  model  to 
introduce  into  the  Lam  model  the  effect  of  ionization  of  r.  eu  t  r  a ; 
gas  molecules  by  current  returning  to  the  probe.  Figure  2.1 
i  1  1  u  s  t  r  a  t  es  sc  h  erria  1 1  cal  1  V  the  ph  vs  i  c  a  1  pr  oc  esses  in  the  sh  e  a  t  h 
region.  The  case  of  a  spherically  symmetric  electron  beam 

emitter  was  considered  and  it  was  found  that  ionization  serv-ec  to 
reduce  probe  charging.  This  was  significant  in  that  experi¬ 
mental  results  then  available  indicated  low  probe  potentials  at 
high  beam  currents .  However  ,  the  study  by  Leadon,  et.al  .  was 
not  entirely  conclusive  in  that  the  PLAS3  program  developed  to 
implement  the  model  failed  to  converge  at  high  beam  currents. 

The  work  on  this  project  was  initiated  with  the  aim  of  deter¬ 
mining  the  origin  of  the  convergence  problems  and  extending  the 
model  to  higher  currents. 

The  further  development  of  the  PL.AS3  code  is  a  step  towards  a 
synthesis  of  spacecraft  charging  models.  FLA S3  is  an 

intermediate  level  program  which  treats  the  basic  physical 
processes  without  dealing  with  the  complications  of  spacecraft 
structural  detail .  Thus  it  can  be  used  to  generate  benchmark 
cases,  for  comparison  with  much  slower  particle  pushing  methods , 
or  with  complex  existing  codes  such  as  MASCAP  and  POLAP  that  do 
not  address  the  low  altitude  ionization  problem. 
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Figure  2.1  Schematic  of  Ionization  Process 


2.  !  The?  PLwS3  Model 


The  mode  1  1  s  essenti  al  1  y  the  same  as-  that  used  in  the  previous 

i  t,4j 

studies.  The  spherical  probe  is  assumed  to  be  surrounded 

bv  a  sheath  of  radius  R  .  At  the  sheath  boundary  the  i  i  <=•  i  d  and 

o 

potential  are  set  to  zero,  and  electron  density  is  equal  to  the 
ambient  plasma  density  n_.  Ambient  ions,  however  .  are  assumed 
to  be  repelled  -from  the  probe  so  that  their  density  is  zero  both 
inside  the  sheath  and  at  the  boundary.  The  sneath  boundary  is 
determined  by  the  condition  that  the*  inward  thermal  electron  flu 
is  equal  to  the  beam  current  I.  ,  thus 

h 


4  H  q  H  r.  1  h  r,e  =  1  b  >.  1  • 

Once  inside  the  sheath,  the  electrons  are  assumed  to  travel  by 
tree  tall  toward  the  probe.  Conservation  ot  energy  then  gives 

J*  m  V*  =  h  m  +  q  *  ,  o. 

On  their  wav  to  the  probe  surface,  electrons  s  trite  neutrals,  at 
a  den s  i  t  v  nr, ,  and  i  on  l  za  1 1  on  oc  c u r  s  w i  t  h  f  r  equ en v 


Where  is  a  velocity  dependent  probability  curve,  normalized 

to  unity,  and  0  is  the  maximum  cross  section  for  ionization, 
around  Z'x  id-'**  cirr-  for  Cv  targets,  with  a  peat  around  ?b  e1  ' . 

The  ion  pairs  thus  created  are  accelerated  away  b-*  the  fie'd  and 
give  rise  to  a  secondary  electron  and  ion  flu*:.  The  number  of 
secondary  pairs  produced  per  second  throughout  the  sheath  is 


where  R j-,  is  the  probe  radius. 
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<3)  A  treatment  equivalent  to  ( 2>  -for  secondary  ions  gives 


r  ^ 


n  s<  r  > 


(&)  >2  dP _ 

(  yj2  +  2g  ($<p) -<tKr))}  ^ 
JR0'  *  mi  J 


1  6  c  > 


Here  the  ions  are  assumed  to  be  swept  out  by  the  motion  of  the 
probe  at  velocity  Vs  so  that  energy  conservation  appears  as 


>4  m  K>2  =  J*  m  vf  +  2a  !♦-*£,  I 

mi 

where  *5  is  the  potential  at  which  the  ion  was  created. 


When  Eqns.  <6 a  -  6c~>  are  substituted  tor  the  RHS  of  (5)  and  the 
boundary  conditions 


•<R0)  »  0 ,  ♦" <R0>  =  0 

are  imposed,  we  arrive  at  a  closed  1 n t egro-di f f er en t 1 al  equation 
•for  ♦,  which  we  solve  by  the  following  numerical  scheme. 


In  order  to  solve  Eq.  5,  we  define 

E  =  -V*  <7> 

and  consider  a  corresponding  difference  equation.  The  sheath 
region  is  divided  into  N  shells  of  width 

Or  =  <r0-rb)/N 

Each  shell  begins  at  a  radius  Pt  =  rb  +  (i-l)4r 
and  has  as  its  midpoint  the  radius  rt  =  rb  ♦  <i-J4>dr 
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we  define 


f 


At  eich  pj  we  define  a.  field  variable  E ,  and  at  each  r  , 
a.  potential  'ti  5  and  a  space  charge  n  j  .  The  potential  is  related 
to  the  field  by 

*1  +  1  '  *1  =  <Ei*2EiM*E,+2>  Jr/4  1  8) 

The  densities  are  found  bv  replacing  the  integrals  in  Eg.  6  bv 
trapezoidal  rule  appro;-  ima  tes  . 

To  include  the  effects  of  ionization  in  cell  i  to  the  charge 
density  in  that  same  cell ,  we  divide  the  cell  into  general 
subcel  is  and  approximate  the  field  as  a  constant  c-  er  each 
r  eg l on .  Th e  int egr  a  1 s  l n  Ea .  6  c  an  than  be  don e  an  a  1 v tic?!  1 y . 

W  i  t  h  the  K  n  ow  ledge  that 


eM+  1 

=  0  , 

'  *  N  +  1  =  0 

We  ar-ri'- 

•e  at  a 

set 

of  M  non-linear 

equ ations  in  N  u n K n own 

fie! 

i  d 

po i n t s . 

These 

we 

solve  by  Newton 

-  R  a.  p  h  s  o  n  iter  a.  t  i  o  n  a  s  f 

O  1  1  O'. 

'..IS 

C'ef  i  ne 

VEi 

,Ey 

'  '  ‘  in 

i  +  l 

-  p‘Ej  -  ezlrrij  , 

i  = )  ,N 

? 

Gu  es  s  an  i  n  i  1 1  a  I  vec  t  or 

»;  o  i  i  i-i  i  i  o  ■* 

'El  ,  E2  . EN  > 

and  modify  the  initial  vector  in  a  direction  which  will  minimize 
tne  residuals  of  Eq.  . 
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Method  of  Approach  and  Result?. 


After  i nves t i ga  t i n  q  both  the  model  and  the  previous  numer  i:al 
implementation,  the  code  was  rewr 1 tten  using  a  different  a  no  m or 
straightforward  method  < Newton -Raph son  •  to  aol'/e  the  equation  a. 
This  program  was  thoroughly  cross-checl ed  with  the  original  for 
cases  where  the  original  converged.  Using  the  new  algor  if  hit', 
we  were  able  to  obtain  solutions  for  arbitrarily  high  beam 
currents.  The  nature  of  the  potential  contours  was  rather 
unusual  and  SUATEK  graphics  capabilities  were  added  to  the 
program  to  facilitate  understanding  of  the  under  lying  pnysics. 
Plots  of  potential ,  electric  field  and  total  and  partial  charge 
distributions  such  as  are  shown  in  Figures  2.2  and  2 . T  were  made 
available.  a  series  of  runs  have  been  made  to  calculate  the 

dependence  of  the  probe  charging  on  the  mode)  parameters,  namei-. 
the  electron  temperature  and  density,  the  mean  ionization  length 
and  ionization  profile,  and  the  beam  current. 

The  effect  of  mean  ionization  distance  X  on  the  potential  of  a 
spac  ec  r af  t  emitting  a  beam  i s  s  h own  in  Figure  2.4.  In  the 
region  where  ionization  is  negligible  1 X  =  Id®  approximately' , 
the  potential  curves  are  flat,  and  the  potential  increases  with 
beam  current.  As  the  ionization  distance  decreases  • i.e. 
ionization  is  more  probable)  the  potential  curve  turns  down,  and 
stavs  constant  approximately  as  the  ionization  distance  X 
decreases  further. 

The  main  result  from  the  parameter  studies  was  that  the  curve  of 
the  probe  potential  versus  current  was  very  similar  to  that 
measured  exper imental 1 y.  For  low  electron  beam  currents,  the 
potential  increases  monoton i cal  1 y  with  the  current.  Beyond  a 
critical  current,  the  returning  electrons  from  fhe  environment 
cause  ionization,  and  the  potential  drops  to  a  lower  ■  alue.1 
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Figure  2.3  Potential  and  Field  Profiles  for  a  Current  Sweep 
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MEAN  IONIZATION  LENGTH  (m) 


E x  t  en  s  1  on  s  of  t  h e  Mode  1 


The  mode!  was  applied  to  two  further  systems.  first ,  the  case 
of  ar<  infinitely  long,  cyl  i  ndr  i  c  al  1  y  symmetric  electron  beam 
emitter  was  considered.  The  necessity  that  the  cyl  i rider  be 
infinitely  long  arises  from  the  fact  that  this  renders  the  mode'1 
one-dimensional  ,  whereas  use  of  further  dimensions  introduces' 
or  ob  1  em  c  omp  1  j  cat  i  on  s  wh  i  c  h  sh  ou  Id  initial  1  y  be  svo  i  oed  . 

Results  from  this  study  were  quite  similar  to  the  spherical  case. 
This  shows  that  the  spherical  nature  of  the  probe  is  not  essen¬ 
tia]  in  determining  charging  characteristics. 

h  case  of  more  immediate  interest  was  the  ion  emitter  Here 

we  needed  to  include  one  further  effect  of  secondary  electron 
emission  from  the  probe  surface  by  primary  ion  collision.  We 
found  that,  for  ion  guns,  ionization  was  less  important  in 
determining  charge  than  for  electron  guns.  In  fact,  the  major 
determiner  was  the  enhancement  of  the  return  current,  mostly  from 
secondary  emission,  to  increase  the  effective  beam  current.  The 
result  of  this  was  to  bring  down  the  magnitude  of  probe  charge, 
as  with  electrons,  but  it  became  clear  that  no  maximum  in  the  1 -V 
curve  would  be  encountered  for  ions. 

nn  analytical  model  was  constructed,  based  on  the  Lam  solution 
and  including  secondary  emission  but  not  sheath  ionization. 

It  was  found  that  this  model  reproduced  the  results  of  the  full 
mode!  quite  nicely  for  environments  of  greatest  interest  to 
experimenters.  This  indicates  that  ionization  may  not  be 
important  for  typical  ion  gun  experiments. 
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2.4  Flight;  Si  mu  I  a  1 1  on 

Studies  have  been  carried  out  to  give  an  indication  of  the 
predictive  usefulness  of  the  model  .  The  goal  was  to 

duplicate  as  nearly  as  possible  the  environmental  conditions  to 
which  an  upcoming  test  flight,  BERT-1,  would  be  subjected. 
Simulations  were  carried  out  at  altitudes  between  180  and  250  km 
for  the  two  electron  gun  models,  spherical  and  cylindrical ,  as 
welt  as  the  spherical  ion  emitter.  The  major  findings  were  the 
following,  assuming  a  maximum  beam  current  of  28  mA: 

a'*  for  electron  guns  the  potent  i  al -current  curve  will  ex  i  bit 
a  local  maximum  for  low  altitudes  but  remain  monotonic  at 
higher  altitudes,  and 

b>  for  ion  guns ,  the  potential  will  reach  the  limiting  value 
of  the  energy  of  the  ion  beam  itself  for  all  altitudes. 

Figure  2.5  shows  the  computed  results  of  spacecraft  potential  as 
a  function  of  electron  beam  current  at  various  altitudes.  The 
non-monoton i c  behavior  of  potential  current  curves  is  apparent. 
At  low  currents,  the  potential  increases  with  beam  current. 

When  foe  current  increases  further,  ionization  occurs  inside  the 
sheath,  and  the  potential  then  turns  around. 

T  n  add i t i on  to  t h e  a tmosph er l c  si mu  1  a t i on ,  a  s  t udy  was  c ar r  l ed 
out  for  electron  guns  at  conditions  similar  to  those  in  a  recent 
ground  based  test  chamber.  In  this  case,  it  was  shown  that  the 
peak:  voltsqe  and  peak  current  could  be  matched  with  reasonable 
choices  of  plasms  density  and  temperature  in  the  model . 
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Figure  2.5  Flight  Simulation  Charging  Characteristics 
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3 . 0  Particle  Trackinci  Studies 

3.1  Particle  Tracking  near  Differentially  Charged  Surfaces 

An  important  aspect  of  the  beam  emitter  problem  is  the  dynamics 
of  the  beam  particles  themselves.  If  a  differentially  charged 
spacecraft  is  placed  in  an  environment  such  that  little  current 
is  drawn  from  the  surrounding  plasma  and  no  sheath  develops, 
Laplace's  equation  holds  and  we  can  calculate  the  potential  at 
all  points  outside  the  spacecraft  knowing  only  the  charge  on  the 
surface.  Once  this  is  known,  calculation  of  the  trajectories  of 
emitted  charged  particles  will  answer  the  question  of  what  frac¬ 
tion  of  an  emitted  beam  will  escape  the  probe,  or  return  and 
thereby  neutralize  some  of  the  charge. 

A  number  of  investigations  were  conducted  to  calculate  both  the 
potentials  around  differentially  charged  spacecraft  and  the 
trajectories  of  emitted  beams.  The  major  results  of  these 
studies  have  been  the  presentation  of  potential  and  particle 
trajectory  plots,  and  the  determination  of  escape  energy  as  a 
function  of  beam  angle  for  several  different  charging  configu¬ 
rations.  Contour  plots  were  used  to  aid  in  the  presentation  of 
these  data.  The  final  objective  of  this  work  is  to  determine 
the  time  dependence  of  sel f -neu tr al i z a t i on  in  a  probe  which  traps 
some  of  its  emitted  beam.  Configurations  studied  included 
quadrupoles,  circular  concentric  disks,  and  spheres  with  a 
circular  cap  —  all  of  which  were  treated  essentially  as  two- 
dimensional  problems. 
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3.1.1  Quadrupole  Problem 


The  potential  tor  two  like  charges  at  a  distance  ' d '  on  either 
side  of  a  central  charge  'i'  is  given  b> 

4><  x  ,y/  =  - f -  -  - ^ ^ - 

<x2  +  y2)^  <<x+d>2  +  y2)^  (ix-d>2  +  y2>  ** 

The  equations  ot  particle  motion  are  given  by 

2 

d  x  =  _  £* 

dt"  &x 


d  y  _  _  i  0 

d7*  iy 

and  were  solved  computationally  using  Hamming's  modified 
predictor-corrector  method.  (subroutine  HPC-G  from  the  IBM 
Scientific  Subroutine  Library) 

Particle  trajectories  were  plotted  using  SUATEK  with  appropriate 
scalinci  to  overlay  on  transparencies  of  the  associated  potential 
contour  plots.  These  contour  plots  were  generated  using  routine 
CONREC  from  the  NCAR  graphics  package.  Figure  3.1  shows  a 
typical  plot  with  four  different  trajectories  for  different  beam 
energies.  Other  parameters  that  were  studied  were  quadrupole 
separation  and  net  charge,  beam  polarity,  and  beam  origination 
and  angle. 
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'z.i.2  Circular  Double  Disk  Problem 

The  potential  tor  a  circular  double  disk  has  been  derived 
analytically^'  as  -follows: 


V(fV)  -- 

+  hi  ('  s  -  - 


/1*MLYTXCAL  CH6CKS  Fo«  Bou*t>ARY  CASES 

V  (■>,*)  =  ifxtay'(i)  *  0-r)[*/(zK«.*//l] 

'**«*)' ,'lj] 


A.W<L  f  >  I 


+  ()- Y)  tao  '[(l~» <-*")  2 


1+  beam  orientation  is  initially  in  a  plane  normal  to  the  disk 
surfaces,  we  can  take  x  =  P  (radius)  and  y  =  z  (normal  to  disk) 
and  consider  the  problem  in  rectangular  coordinates. 
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The  inner  disk  to  outer  disk  radius  ratio  a,  and  the  potential 
ratio  Y  provide  parameters  that  are  similar  to  the  quadrupole 
problem.  The  main  complexity  is  due  to  the  integral  which  must 
be  computed  for  the  potential  at  each  location  P.2. 

Figure  3.2  shows  overlaid  potential  and  trajectory  plots  for 
positively  charged  inner  and  outer  disks,  with  6.1  radius  ratio. 
Ion  beams  with  three  different  energies,  originating  at  a 
distance  of  6.2  above  the  disks  are  considered.  The  saddle 
point  introduces  the  question  of  beam  trapping  for  many  of  the 
studies  described  in  this  section. 


3.1.3  Circular  Cap  Problem 


The  problem  is  to  track  particles  in  the  vicinity  of  a  spherical 
probe  with  a  differentially  charged  surface. 


\ 

\ 

v 


' v  =  y  volts 


1  or  0  vol ts 


The  potential  is  calculated  by  noting  that  the  series 

oo  n+ 1 

LKr,6>  =  Z  (£)  Of,  Pn(cos6)  <  1) 

n=0V.r/ 


solves  Laplace's  equation,  where  a  is  the  probe  radius. 
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CIRCULAR  O0UBLE  DISK 


The  coefficients  are  easily  found  since 


and 


V<a,©> 


{ 


r 

1  (or  0) 


$  <  P 

e  >  9 


«n 


<n  +  S> 


ft 

f  <*» 
o 


Pn(cos$)  d(cos(» 


by  orthogonality. 


Program  TRKCOEF  calculates  these  coefficients.  It  is  run  inter¬ 
actively  and  requires  the  following  inputs: 


NTERM  -  #  of  terms  in  the  series 


RRR  =  cutoff  distance  for  "nearness"  to  the  probe. 

(good  convergence  is  obtained  for  RRR  =  1.1  when  NTERM=38> 

B  =  the  angle  £  in  degrees 

G  —  the  charge  on  the  cap 

The  integrals  involved  are  evaluated  analytically;  the  limit  on 
NTERM  is  around  50  due  to  the  difficulty  of  calculating  high 
order  Legendres. 

Whether  the  charge  on  the  lower  portion  of  the  probe  is  1  or  0  is 
determined  by  the  parameter  GAB  in  subroutine  CHECK.  Program 
TRKPLOT  sets  up  a  file  for  contour  plotting.  TRKCOEF  saves  the 
coefficients  on  TAPE5  and  must  be  run  (or  TAPE5  saved)  before 
other  programs  are  run. 
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Program  TRK10NS  uses-  HPCG  to  track  particles  from  a  previously 
generated  potential.  It  allows  various  starting  positions,  on 
or  off  the  probe,  and  various  beam  energies  and  angles.  They 
are  defined  as  follows: 


*  start  distance  =  0  wi 1 1  start  the  particle  at  r  =  RRR. 

Every  timesteps,  the  program  asks  whether  you  wish  to  stop. 

Thus,  in  its  present  form  it  has  to  be  run  interactively. 

The  status  of  this  problem  was  that  the  “hot  spot"  did  not  seem 
to  trap  the  beam  when  it  is  released  from  a  point  near  the  craft. 
This  means,  at  least,  that  the  situation  is  more  complicated  than 
was  originally  imagined. 
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3.2  Evaluation  of  Geomagnetic  Field  Models 

I n  tr  oduc  t  i  on 

An  important  -feature  in  the  application  of  particle  beam  emitters 
is  the  accuracy  with  which  the  beam  can  be  delivered  from  the 
spacecraft  to  the  selected  target.  Since  the  beam  is  charged, 
its  trajectory  is  strongly  influenced  by  the  geomagnetic  field  it 
must  traverse.  Thus  accuracy  of  timing  is  contingent  on  an 
accurate  knowledge  of  the  geomagnetic  field  itself.  This  study 
determines  directly  the  effect  of  uncer t an t i es  in  the  geomagnetic 
field  model  on  targeting  by  calculating  the  trajectories  of  beam 
particles  with  various  geofield  models.  Since  the  discrepency 
within  different  models  is  expected  to  be  at  least  as  large  as 
the  difference  between  the  “best"  model  and  the  true  field,  the 
study  should  provide  a  rough  upper  bound  on  the  accuracy  with 
which  beams  can  presently  be  targeted. 


Model s 

The  models  are  based  upon  fits  of  experimental  data  to  an 
expansion  of  spherical  harmonics.  The  coefficients  are  referred 
to  a  specific  year  and  correction  terms  are  sometimes  provided  or 
subsequent  years.  The  study  was  conducted  for  the  year  1980. 

In  addition  to  yearly  variation,  there  is  also  a  dependence  in 
the  field  on  time  of  day  and  sunspot  activity.  These  effects 
were  investigated  by  using  additional  sets  of  coefficients, 
consisting  of  correction  terms  to  the  reference  I6RF80  set,  for 
0900  and  1300  hours  and  for  noon  at  sunspot  minimum.  The  models 
used  and  the  size  of  the  coefficient  sets  are  given  in  Table  3.1. 
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Table  3.1  Geomagnetic  Field  Models 


Model 

1) 

IGRF80 

2) 

Barrac 1 ough 

80 

3) 

Barrac 1 ough 

75 

4) 

Magsat 

5) 

WC80 

6> 

Sq  mods  <  to 

IGRF80) 

Order  <NMAX  =  ORDER  +  1>  j 


10X  10 
10X  10 
12x  12 
13x  13 
12x  12 


Simulations  were  carried  out  tor  two  particle  types,  two  starting 
altitudes,  three  initial  beam  energies  and  three  aiming  direc¬ 
tions.  Thus  r,  total  ot  36  trajectories  were  calculated  tor  each 
model.  Trajectories  were  Followed  through  one  gyro  rotation  and 
stopped  20  km  below  the  starting  altitude.  The  I6RF80  trajec¬ 
tory  was  used  as  the  target  detinition,  and  deviation  ot  other 
trajectories  From  this  one  was  calculated  as  the  nearest  approach 
ot  a  target  and  test  trajectory,  without  regard  to  correlation  in 
time. 


Procedure 

In  order  to  add  a  set  ot  coetticients  to  the  existing  system  ot 
programs,  the  Following  sequence  ot  runs  is  carried  out: 

1>  The  program  CNGMAG  is  used  to  convert  the  Gauss  normalized 
coetticients  to  Schmidt  normalized  Form  and  to  correct  the 
coetticients  to  the  desired  epoch.  The  correction  terms  are 
input  along  with  the  coetticients  and  currently  extend  to  second 
order  in  time.  The  program  creates  a  data  block  containing  the 
renormalized  coetticents  to  be  injected  into  subroutine  NEWMG10 
ot  the  trajectory  generating  program. 
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2)  The  simulator^)  was  provided  by  D.  Smart  <AF6L;>  and  was 
transferred  to  a  program  library  called  T.JPL.  There  are  several 
copies  therein  of  the  field  generating  subroutine  NEWMG18.  The 
highest  order  spherical  harmonic  to  date  is  N=14.  If  a  higner 
order  is  desired,  the  program  CDNMA6  will  write  a  new  NEWMG10 
subroutine.  Some  problems  in  the  format  statements  of  CDNMaG, 
causing  intermediate  zeros  in  constants  to  be  eliminated,  have 
been  encountered  and  care  should  be  exercised  in  any  future  use. 
If  the  dimension  of  the  new  coefficient  set  is  14  or  less,  one  of 
the  existing  NEWMG10's  in  TJPL  can  be  used  directly. 

3)  The  data  statements  written  by  CNGMAG  are  injected  into 
NEWMG18  and  the  program  library  is  updated. 

4)  The  simulation  is  carried  out  for  a  given  starting 
configuration  and  field  model  using  TJPL.  Figure  3.3  shows  a 
typical  trajectory  for  a  few  gyro  rotations. 

5)  Results  from  this  and  previous  IGRF88  calculations  are 
compared  using  program  TJCALC,  which  finds  the  minimum  distance 
discrepency  between  test  and  reference  trajectory. 

Resu l ts 

The  longest  distances  travelled  occurred  for  high  energy  lithium 
ions  at  high  altitude  and  were  on  the  order  of  780  km.  The 
shortest  were  for  low  energy  protons  at  low  altitude  and  were 
around  15  km.  The  miss  distances  from  the  IGRF  80  trajectory, 
as  a  percent  of  total  distance  travelled,  were  0.887.  for  WC  38 
and  0.817  for  MGST  4181  at  high  altitudes,  and  0.147.  for  WC  88 
and  0.0257.  for  MGST  4181  at  low  altitude.  The  daily  solar 
corrections  added  to  the  IGRF  80  model  gave  target  deviations  or 
the  same  order  of  magnitude. 
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4 .  ®  Atmospheric  Density  and  Winds  -from  Accelerometer  Data 

Triaxia)  accelerometers  capable  of  measuring  accelerations  in 
the  tenths  of  micro-g's  are  released  in  rocket  launched  spinning 
spheres.  Once  the  sensor  system  is  not  subject  to  the  large 
propulsive  forces,  the  sensitive  devices  are  uncaged  and,  being 
in  continuous  free  fall  in  the  earth's  gravitational  field,  are 
able  to  detect  the  remaining  forces  acting  on  the  sensors. 

These  forces  are  primarily  due  to  the  drag  of  the  sphere  in  the 
earth's  atmosphere,  but  may  include  centrifugal  forces  arising 
out  of  module  spin  and  precession.  Since  the  objective  of  these 
experiments  is  the  study  of  atmospheric  structure  in  the  lower 
thermosphere,  eccentric  forces  are  an  undesirable  contaminant, 
and  are  eliminated  by  design  by  locating  the  sensors  close  to  the 
dynamic  center  of  the  spinning  spheres.  Any  residual  imbalanced 
signals  are  attenuated  by  special  notch  filters  when  the  expected 
contaminant  frequencies  are  known  e.g.  spin  frequency  for  spin 
axis  2-sensors.  Other  undesirable  modulation  signals  must  be 
totally  removed  by  computational  filtering  prior  to  any  detailed 
processing  for  analyses  of  atmospheric  density  and  wind 
structure. 

In  the  case  of  the  rocket  launched  experiments,  the  sensor  unit 
resides  within  a  sphere  that  is  released  after  launch,  and 
includes  the  telemetry  for  transmitting  the  sensor  signals.  The 
trajectory  is  geographically  "local",  and  analyses  have  been 
performed  in  launcher  referenced  coordinates  in  which  the  spin 
axis  orientation  stays  inertial 1y  fixed.  Earth  curvature  must 
be  accounted  for  when  equating  upleg  and  downleg  densities  and 
wind  vectors,  and  can  be  accomplished  by  using  actual  trajectory 
altitudes  at  all  times.  Figure  4.1  shows  the  coordinate  systems 
which  are  appropriate  for  rocket  sphere  studies.  The  X"  system 
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is  in  the  spinning  sphere  reference  frame,  with  2“  the  spin  axis 
direction.  The  X'  and  X  systems  are  non-spinninq  and  embedded 
in  the  trajectory  plane,  with  2/<=Z“)  coincident  with  the  sphere 
axis,  tilted  forward  at  angle  a  relative  to  launcher  vertical  2. 


Figure  4.1  Coordinate  Systems  for  the  Rocket  Sphere 


Important  functions  that  the  data  processing  system  performs  are 

1)  Telemetry  data  compaction  and  splicing 

2)  Sensor  calibrations 

3)  Signal  demodulation  and  conditioning 

4)  Trajectory  evaluation  and  integration 

5>  Maintenance  of  an  evolving  Resource  File 
6)  Resolution  of  sensor  orientation  and  transverse  and 
in-track  winds 
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7) 

8) 

9) 

10) 
1  1) 


Resc  1  u  t  i  on 
prof i 1  es 
Tables  and 
batch  mode 
Compar i son 
Eval uat i on 
Studies  of 


o f  consistent  density  and  temperature 

graphical  presentations  in  interactive 

with  other  concurrent  experiments 
of  reliability  of  results 
auxiliary  associated  phenomena 


or 


A  schematic  of  the  program  f 1 ow  and  intermediate  tiles  is  shown 
in  Figure  4.2,  starting  from  the  DECGMM  tape  and  calibration  and 
trajectory  data.  In  addition  to  the  final  density  and  wind 
profiles  in  file  SPIR,  tables,  plots  and  supporting  analyses  are 
end  products.  The  scope  of  each  functional  area  is  now 
discussed.  The  AC13  sphere  experiment  at  Esrange  on  Dec . 1 ,  1980 
serves  as  the  prototype  for  the  sample  results  in  this  report. 


4.1  Telemetry  Data  Compaction  and  Splicing 

The  12-bit  PCM  coded  sensor  signals  are  available  over  the 
duration  of  the  flight,  after  sphere  release.  It  is 
advantageous  if  the  complete  data  are  available  on-line  during 
the  initial  splicing,  calibration,  and  sensitivity  and  phase 
adjustment  studies.  PACKER  routines  available  in  the  SUNY 
program  library  allow  all  the  data  for  four  X  and  Y  channels, 
plus  three  Z  and  a  Nutation  channel  over  the  300+  second  span  of 
the  experiment  to  be  resident  on  about  2000  PRU's  of  perm  file. 
The  DECQMM  tape  is  checked  for  sync  losses,  and  the  packed  file 
is  augmented  from  alternate  telemetry  channels  for  any  missing 
data.  Time  track  errors  of  nearly  1  second  between  channels 
from  different  radars  have  been  identified  during  the  splicing 
process;  thus  such  procedures  have  enhanced  quality  control  over 
the  data  base.  Once  sensor  calibrations  and  signal  demodulation 
results  are  completed  and  saved  on  reduced  files,  this  large 
packed  file  need  not  be  resident  on  permanent  file. 
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Figure  4.2  Schematic  of  the  Rocket  Sphere  Processing  System 
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4.2  Sensor  Calibrations 

Uncaging  of  the  sensors  after  sphere  release  initiates  a  damped 
resonant  response  in  all  the  accelerometer  bimorphs.  These 
resonant  frequencies  are  invaluable  in  fixing  the  calibrations, 
because  the  vibrations  are  produced  within  the  centrifugal  force 
field  being  experienced  by  the  sensors,  and  the  resultant 
stiffness  and  sensitivity  calibration  estimates  improve  on  the 
•values  obtained  during  lab  tests.  Correction  factors  are 
applied  by  the  researcher,  based  on  the  ratio  of  in-flight  and 
ground  based  resonant  frequencies.  Optimized  function  fitting 
techniques  are  used  to  converge  on  a  best  fit  solution  for  each 
damped  oscillation  frequency.  The  Z-axis  cases,  where  the 
resonance  is  superposed  on  the  precession  modulated  signal,  have 
also  been  implemented  and  solved  using  the  appropriate  extended 
expressions.  A  different  type  of  procedure  is  necessary  to 
determine  a  calibration  correction  that  accounts  for  the 
attenuated  slow  time  constant  response  of  the  Z-sensors  as  the 
sphere  experiences  exponentially  higher  densities  on  re-entry. 
Z-signal  time  constants  are  determined  from  log  plots,  and  are 
applied  to  lab  tests  that  were  performed  using  various 
exponential  driving  functions.  A  pre-calibration  package  of 
results  and  supporting  plots  has'e  been  provided  for  each  sphere 
flight. 


4.3  Trajectory  Evaluation  and  Integration 


In  general  the 

the  drag  forces 

bv  P  C  .  Area 
d  ~~ 2m~ 
relative  to  the 


trajectory  velocity  is  the 
experienced  by  the  sphere, 
where  is  the  veloci 

atmosphere,  it  can  be  seen 


greatest  contri 
Since  drag  i 
ty  of  the  spher 

that  inaccurac 


butor  to 
s  given 
e 

i es  in 
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V^  are  amplified.  A  simplified  sphere  trajectory  calculating 
routine  was  installed,  a)  to  create  good  quality  quick  look 
trajectory  velocities  before  refined  results  become  available, 
and  b)  to  evaluate  initial  trajectories  provided  by  the  range. 
Launcher  referenced  trajectory  velocity  components  are  formatted 
and  entered  at  2  to  5  second  intervals  in  the  experiment  Resource 
File.  Velocities  at  any  time  or  altitude  upleg  or  downleg  are 
obtained  by  cubic  interpolation. 

4.4  Evolving  Resource  and  Profile  Files 

In  addition  to  trajectory  information,  the  concept  of  the 
Resource  File  is  to  collect  all  initial  parameters  and 
intermediate  results  that  are  required  for  processing  of 
subsequent  stages  of  the  analysis.  Thus,  sphere  mass,  cross- 
section  area,  original  sensor  calibrations,  etc.  are  initially 
entered.  The  nominal  spin  rate  is  established  early  and  is 
entered  manually.  Calibration  adjustmen ts  to  give  continuity 
and  consistency  are  determined  later,  and  are  then  entered  and 
called  by  subsequent  programs.  Figure  4.3  shows  a  sample 
segment  of  the  Resource  File  for  the  AC13  flight.  This 
presentation  provides  the  investigator  and  the  program  user  with 
an  up-to-date  status  of  the  investigation  for  each  experiment. 
Profile  wind,  density,  and  related  solutions  per  altitude  step 
form  a  much  larger  data  base,  and  are  not  stored  on  the  Resource 
File,  but  on  successively  augmented  binary  files  for  ready 
processing  and  plotting.  Since  many  rocket  sphere  flights 
continue  to  be  studied  through  improved  multi-stages  of  proces¬ 
sing  which  generate  intermediate  profiles,  a  logistics  control 
procedure  must  be  maintained.  Figure  4.4  shows  a  recent  chart 
which,  using  established  input  from  the  Resource  File,  permits 
any  of  the  program  modules  to  be  rerun  with  minimum  effort. 


AC  1  3 


ESRANGE 


12/01/80  540.00  SEC. 

SMASS , SAREA  (KG.  M««2) 

10.262,  0.05199 

ALAT.FLAZ  (LATITUDE,  FLIGHT  AZIM.) 

68.,  -10.5 

FRATE  (FRAME  RATE  PER  SEC.) 

.008 

CALV  (CALIB.  VOLTS  PER  COUNT) 

.009775 

RPS  (SPIN  RATE,  OFFSET.  DRIFT  REVS.  PER  SEC.) 

5.93340,  -.0008,  .20E-4 
ALPHA, BETA  (AXIS  TILTS  OEG. ) 

14.  .  0. 

MXY  ( 1 *USE  X-SENSOR  2*USE  Y-SENSOR) 

2 

TXY  (3-2.  2-1,  1-2,  2-3,  3-4  CROSS-OVER  TIMES  X  AND  Y) 

82. . 94. .285. .312. .331 .5,  80 . . 93 . 5 , 286 . . 31 3. 5 , 333 . 5 

CALXY  (ORIG.  CAL.  M/S*. 2  PER  VOLT  FOR  X1,X2,X3,X4,  Yt,Y2,Y3.Y4) 

6.8375E-6, 1 . 8071 E-4 , 4 . 721 6E-3 , .1451 
9 . 233E-6 , 2 .489E-4 ,6 . 441 OE-3 . . 1926 

CRXY  (EST.  SENS.  RATIOS  X 1 /X2/X3/X4/Y4 .  Y1/Y2/Y3/V4/Y4) 

26.45,25.90,31.16,1.23,  30.25,24.75,29.92,1.00 

CPXY  (EST.  PHASE  OIFFS  X1-X2-X3-X4-X4.  Y1-Y2-Y3-Y 4-X4 ) 

-7. .45. ,-56. .0.0.  -3. .42.2,-60.5,-88.2 

PHC  (PHASE  BIAS  X-SENSOR  AND  Y-SENSOR) 

270. . 270. 

PHS  (PHASE  ADJUSTMENTS  X1.X2.X3.X4.  Y1.Y2.Y3.Y4) 

18. .  1 1 . .56. .0. ,  -66.8,-69.8,-27.5,-88. 

CPH  (LOW  ALT.  SPIN-OOWN  PHASE  ADJ.  »  B* ( A-A LT ) • »2 ) 

85.,  .045 

CAL  (REV.  CAL.  M/S»»2  PER  VOLT  FOR  X1.X2.X3.X4,  Y1.Y2.Y3.Y4) 

7 . 1 90E-6 . I . 902E-4 . 4 . 926E-3 . . 1 535 
8 . 42SE-6 , 2 .550 E-4 .6 . 3 1 0OE-3 , .1888 
PREF  (PRECESSION  FREQUENCY) 

.7248 

TTST  (START  TIMES  FOR  21, Z2, 23, NUT  SENSORS) 

300. . 310. .325. .340. 

TTND  (END  TIMES  FOR  Z1,Z2,Z3,NUT  SENSORS) 

322. . 335. .350. .359.9 

TRANS  (21-22,  22-Z3,  23-NUT  TRANSITION  TIMES) 

318. . 322. .331 . .335. .350. . 350. 

TCROSS  ( Z 1 .  22,  23,  NUT  CROSS-OVER  TIMES) 

320. . 333. .350. .350. 

BG  (BIAS-COUNTS  AND  GA I N-M/S* *  2/VO LT  FOR  Z 1 , Z2 , Z3 , NUT ) 

510. . 5.942E-3.  510. 71 , .05608,  510.36,17012,  512. ,0. 

BIAS  (REV.  BIAS  COUNTS) 

510. . 510.41 .510.22.512. 

2CAL  (REV.  CALS.  M/S«*2/V0LT  FOR  2 1 , Z2 , Z3 , NUT ) 

5 . 838E-3 , 5 . 72 1 E-2 , .6996,0. 

2TST.2TND  (START, END  OF  Z-NUT-SENSOR  PROCESSING) 

300..  345. 

BOT.HNO  (BOTTOM  ALT.  FOR  XYUP,  END  ALT.  FOR  PROCESSING) 

85.,  55. 

ENOROK  •**.*.«.*,.***•»*«••  END  OF  ROCKET  RECORD 


Figure  4.3  Evolving  Resource  File  (sjunple  segment) 
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PROCEDURE 

PROGRAM 

DEMODULATE  Z-SENSOR  SIGNALS  USING  KNOWN 
PRECESSION  FREQ.  TO  GIVE  COUNTS/SENSOR 

ZEEDEM 

AC77TAP 

ZEEDEM 

OETERMINE  OPTIMUM  BIAS  AND  GAIN  RATIOS 
BETWEEN  Z-NUT  CHANNELS,  USING  22  BIAS 
ESTIMATE  FROM  INSPECTION  OF  PLOTS. 

AOUUST  Z-CALS  LEAVING  Z2»Z3»24  PRODUCT 
UNCHANGED.  APPLV  Z-CALS  AND  BIASES  TO 
CREATE  2ACC  FILE  FROM  ZEEDEM  COUNTS 

ZACC 

ZEEDEM 

ZACC 

ESTABLISH  REVISED  CALS  AND  PHASE  ADJUST 

XYAMPH 

AC77TAP 

NONE 

TO  GIVE  CONTINUITY  OF  XYACC  AND  PHASE 

FOLLOWED  BY 

BETWEEN  X-SENSORS,  Y-SENSORS.  AND  X-Y 

xytran 

AC77TAP 

XYTRAN 

COMPLETE  XYACC  RUN  USING  FINAL  GAIN 

AND  PHASE  CALS  WITH  ZERO  PHASE  ADJUST 

XYACC 

AC77TAP 

XYACC 

SMOOTHING  AND  1/4  KM.  INTERVAL  DATA 
GENERATION  FROM  FULL  XYACC  FILE. 

THIS  IS  THE  ZERO  PHASE  ADJUST  FILE 

XYAPUD 

XYACC 

APUD1 

INTERACTIVE  SOLUTION  FOR  RPS.  PHC . 

BETA.  ALPHA.  SPIND,  SALT.  DELRPS 

TO  MATCH  UPLEG-COWNLEG  CROSSWINDS 
(UPDATE  SPIPH1  AFTER  FIXING  ALPHA 

USING  X2ACC,  XZRAW .  AND  ERRAN) 

SPIPH 

APUD1 

SPIPH1 

GENERATE  1/4  KM.  INTERVAL  XZACC  FILE 

FROM  ORIGINAL  ZACC  FILE.  USING  APUD1 

AND  OPTIONAL  SP 1  PHI  FILE 

XZACC 

ZACC 

APUD1 

(SPIPH1) 

XZACC 

USE  SUAMRG  TO  MERGE  XZACC  FILE  WITH 

APUD1  FILE  GIVING  COMMON  DATA  FILE 

SUAMPG 

XZACC 

APUD2 

USE  OVERLAPPING  XY  AND  Z  DATA  TO 

SOLVE  FOR  FURTHER  DOWNLEG  WIND 

PROFILE,  AND  TO  RE-ESTIMATE  ALPHA 

XZRAW 

SPIPH1 

APUD2 

XZRAW 

SMOOTH  UP-DOWN  PHASE  AND  DOWNLEG  WIND 
SOLUTIONS  BEFORE  DENSITY-TEMPERATURE 
CALCULATIONS,  AND  FOR  FINAL  PLOTS 

5UATEK 

XZRAW 

SP I PH2 

ERROR  ESTIMATION  PROGRAM  MAY  BE  RUN 

HERE  TO  OBTAIN  SENSITIVITY  ESTIMATES 

FOR  EFFECTS  OF  CALS,  TRAJ  VELS,  ALPHA, 

AND  PHASE  ON  RHOCD.  VXW,  VYW 

ERRAN 

SP I PH2 
APUD2 

NONE 

ITERATIVE  RUN  TO  SIMULTANEOUSLY  SOLVE 

FOR  RHO,  TEMP,  CO  CREATING  THE  FINAL 
SOLUTION  FILE  FOR  DENSITY.  TEMPERATURE, 

AND  WINDS  (N-S  AND  E-W  USING  FLT.aZ.) 

SPIRT 

SPIPH2 

APUD2 

SPIR 

DETERMINE  RICHARDSON  NUMBER  VS.  ALT. 
(TURBULENCE) 

RICHNO 

SP1H 

RICHNO 

MIND  VECTOR  PROFILE  PRESENTATION 

WIVEC 

SPIR 

WIVEC 

_ 

Figure  4.4  Processing  System  for  Density,  Temperature,  Winds 
(programs  to  calibrate,  plot,  tabulate  not  shown) 


4.5  Signal  Demodulation  and  Conditioning 


The  equatorial  plane  <X-Y>  sensors  are  spin  -frequency  modulated, 
with  the  amplitude  of  the  signal  a  direct  measure  of  the  lateral 
drag  forces  experienced  by  the  sphere.  Sine  wave  function 
fitting  is  employed  to  determine  the  modulating  lateral  force 
amplitude  and  phase.  Calibration  sensitivities  between  sensor 
channels  usually  turn  out  to  be  somewhat  inconsistent,  so  that 
sensed  acceleration  values  between,  say,  X3  and  X4 ,  exhibit  a 
jump  discontinuity.  A  program  package  was  developed  to  adjust 
all  initial  calibration  sensivities  so  that  the  net  percent 
change  for  all  sensors  is  minimized,  while  imposing  continuity 
in  the  transition  regions,  and  overlap  between  X  and  Y  sensors. 
Care  is  needed  in  this  trimming  operation,  to  avoid  interpreting 
transition  discontinuities  as  density  profile  features. 

The  Z-axis  sensors  are  not  sensitive  to  the  sphere  spin,  but  the 
signals  are  modulated  at  the  precession  frequency,  which  here 
must  be  rejected  and  not  demodulated.  The  resulting  signal  for 
each  channel  starts  from  near  zero  and  builds  up  exponentially, 
with  the  channel  gains  set  to  overlap  while  providing  maximum 
dynamic  range,  as  for  the  X-Y  sensors  above.  Gain  optimization 
is  again  employed  for  continuity,  but  special  attention  is 
required  to  estimate  and  subtract  a  built-in  non-zero  bias. 
Figure  4.5  is  a  composite  plot  of  the  final  clean  accelerations 
sensed  by  the  Z-sensor  down l eg,  and  the  X-Y  sensors  upleg  and 
downleg.  For  the  X-Y  sensors  in  particular,  the  high  frequency 
spin  modulated  acceleration  signal  6Hz>  results  in  excellent 
bias  free  data.  Demodulation  is  performed  at  ^  spin  cycle,  or 
at  less  than  0.1  sec  steps.  At  the  lower  altitudes  where  the 
sphere  fall  velocity  exceeds  1  km/sec,  this  still  allows  solu¬ 
tions  about  100  m  apart.  Analyses  to  date  have  been  conducted 
conservatively  at  1/4  km  intervals,  but  this  extension  to  higher 
resolution  can  be  readily  implemented,  for  more  detailed  studies 
of  atmospheric  wave  and  structure  patterns. 
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Figure  4.5  Z-sensor  and  XY-sensor  Accelerations 


4.6  Sensor  Orientation  and  Winds 

Use  of  a  nominal  sphere  spin  rate  results  in  a  stable  solution 
for  acceleration  as  a  function  of  altitude.  The  phase  or 
apparent  direction  of  the  acceleration  force  on  the  spinning 
sphere  is,  however,  extremely  sensitive  to  any  offset  between 
true  and  nominal  spin  rate.  The  difference  is  evident  as  a 
runaway  condition  where  the  apparent  peak  of  the  acceleration 
does  not  remain  locked  in  the  general  direction  of  the  trajec¬ 
tory.  An  effective  procedure  has  been  devised  to  optimize  for 
true  spin  rate  and  reference  phase,  so  as  to  result  in  transverse 
wind  velocity  solutions  that  are  as  close  as  possible  on  upleg 
and  downleg.  The  transverse  wind  is  given  by 

V  =  1  EACC  sine  <4. 1) 

w  rovT 

where 

R< h>  =  V_  p  C.  Area  (4.2) 

R  d  2m 

EACC  is  the  X-Y  or  equatorial  sensed  acceleration, 

6  is  the  phase  to  be  optimized, 

Op  is  roughly  the  trajectory  velocity,  and 

P  approximates  the  true  density  profile  after  a  few  iterations. 

The  phase  e  obtained  originally  using  the  nominal  spin  rate 
becomes  modified  according  to 

S'  =  e  -  <T  +■  ♦  (4.3) 

where  «  is  the  spin  rate  vernier  offset,  T  is  trajectory  time,  ♦ 
i=  the  reference  phase  adjustment,  and  0'  is  the  corrected  phase. 

Upleg  and  downleg  transverse  wind  sensitivities 

iVvw  and  Wvw 


are  then  used  to  adjust  <  and  $  for  closest  upl eg-downl eg  match. 
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Low  altitude  spin  down  is  also  estimated  in  this  manner  by 
optimizing  an  altitude  dependent  polynomial,  and  in  -fact 
introduces  the  possibility  of  investigating  spin  deceleration 
models.  A  significant  achievement  was  the  identification  and 
quantification  of  an  internal  damping  spin  deceleration  for  AC13, 
showing  up  as  a  parabolic  phase  deviation. 

Figure  4.6  shows  the  development  steps  at  this  stage.  Figure 
4.6(a)  shows  typical  upl eg-down 1  eg  phase  solutions  for  a  good 
fixed  nominal  spin  rate.  Figure  4.6(b)  shows  the  final  revised 
spin-phase  after  the  optimization  process,  which  gives  the 
crosswind  match  of  Figure  4.6(c). 

A  crucial  part  of  the  atmospheric  profile  analysis  is  the 
determination  of  the  forward  tilt  a  of  the  sphere  axis  in  the 
trajectory  plane.  The  basic  sensed  acceleration  equations  in 
the  non-rotating  coordinate  system  of  the  sphere  fixed  in  the 
trajectory  plane  are 

XYZ  =  [tf.  ♦  3  E?(P)  £(oc)  P  C  .  Areal  (4.4) 

acc  t  w  R  d  — — — 

2m 

where  subscripts  t  and  w  refer  to  the  trajectory  and  wind 
respectively,  A(a)  ,  BO)  are  spherical  coordinate  transformations 
corresponding  to  sphere  axis  tilt  in-,  across-  the  trajectory 
pi ane . 

The  cross  tilt  £  is  almost  indistinguishable  from  a  small 
difference  in  spin  rate,  in  its  effect  on  the  sensed  accel¬ 
erations.  When  sun  sensors  are  installed  on  daytime  flights, 
this  parameter  can  be  further  studied  along  with  spin-down. 

The  procedure  for  establishing  a  consists  of  assuming  low  winds 
in  the  60-80  km  region  where  X-Y  as  well  as  Z  sensor  accelera¬ 
tions  are  available,  and  then  applying  Eq.  4.4.  Once  ct  is 
fixed,  the  above  equation  is  also  applied  as  follows: 
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a>  In  the  region  above  80  km,  in-track  winds  and  density  are 
solved  simultaneously  using  upl eg-down 1  eg  X-Y  data. 
b>  In  the  88  km  and  lower  region,  in-track  winds  are  solved 
using  downleg  X  and  2  data  tor  the  -fixed  a,  and  then  with 
winds  known,  density  is  solved  using  downleg  X  data  alone. 

Figure  4.7  shows  the  final  estimated  in-track  wind  profile. 
Problems  of  reliability  of  2-sensor  calibration  and  filter  time 
response,  and  therefore  uncertainty  in  the  80  km  transition 
region  have  been  encountered.  Overlapping  upl eg  X-Y  and  downleg 
Z  data  are  usually  available  in  the  80-160  km  region,  but 
unfortunately  the  2  data  are  least  reliable  here.  Improved 
recovery  of  the  true  2-accel erat i ons  should  be  feasible  by  a 
procedure  that  inverts  the  sensor  response  transfer  function. 


4.7  Resolution  of  Density-Temperature  Profiles 

Equation  4.4  shows  that  once  trajectory,  wind  velocity  vectors, 

and  the  sphere  orientation  are  known,  PCd  is  also  solved 

simultaneously  and  forms  the  next  stage  irreducible  data  base. 

C,  is  a  function  of  Reynolds  number<R  >  and  Mach  number <M> ,  which 
a  e 

in  turn  depend  on  temperature  and  the  density  profile.  6iven 
PC^  and  profiles,  this  interrelationship  necessitates  an 
iterative  solution  procedure  as  shown  in  Figure  4.8. 

The  temperature  profile  Tm  is  obtained  by  integrating  along  the 

density  profile  under  the  assumptions  of  the  ideal  gas  law  and 

hydrostatic  equilibrium.  The  drag  coefficient  C.  is  obtained 

(6)  ^ 

using  Whitfield's  model  for  1 ow  R  ,  and  provides  smooth 

e 

transition  between  the  free  molecular  flow  and  the  continuum 
regions/ Figure  4.9  profiles  the  density  ratio  tr.  the  US 
Standard  Atmosphere  <1976),  and  the  associated  temperature,  that 
was  obtained  for  AC  13. 
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0.0  PC  13  8S.  .045  M=5. 932600  PHC=270.0 
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Figure  4.7  Composite  Estimate  for  In-track  Wind  Profile 


Figure  4.8  Density-Temperature  Profiles  from  f and  Vr 

(the  nonlinear  dependence  of  C<j  on  these  parameters 
necessitates  an  iterative  solution  procedure) 


Solved  Density  Ratio  (to  USSA76),  and  Temperature  Profiles  (AC13) 
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4.8  Tabular  and  Graphical  Presentations 

The  processing  system  propagates  augmented  profile  solution 
results  in  SUATEK  binary  format,  and  tabulations  of  desired 
results  are  accomplished  by  simple  print-out  programs.  The 
Tektronix  plotting  capabi 1 i ty  inherent  with  the  files  permits 
intensive  interactive  operation  and  evaluation  of  results  of 
successive  processing  programs.  This  capability  has  permitted 
the  incorporation  of  many  of  the  analytical  features  described 
before,  which  otherwise  would  be  extremely  tedious  in  a  batch 
processing  mode.  The  present  system  allows  for  repeated 
refining  passes  after  the  investigator  evaluates  results,  e.g. 
revised  axis  tilt,  calibration  modification,  spin -down  model 
adjustment,  consideration  of  improved  drag  coefficient  model,  or 
installation  of  other  processing  modules. 

We  include  here  a  tabular  summary  <Table  4.1)  of  atmosphere  and 
wind  profiles.  The  latter  is  also  depicted  in  a  specialized 
plot  (Figure  4.10)  which  combines  altitude  and  wind  vector  data. 


4.9  Comparison  with  Other  Concurrent  Experiments 

Robin  sphere  and  chemical  release  experiment  results  become 
available  as  the  final  results  for  all  aspects  of  a  test  or 
campaign  are  collected  for  coordination.  Good  density  and  wind 
profile  estimates  are  usually  provided  from  Robin  sphere  analyses 
performed  around  60-70  km.  altitude,  which  is  where  the  rocket 
sphere  spin  deceleration  and  therefore  the  wind  and  density 
profile  is  least  reliable.  Improvement  of  the  spin  deceleration 
model,  and  even  the  sphere  axis  tilt  may  be  warranted  if  the 
consensus  from  the  other  tests  so  dictates.  These  test  wind 
velocities  are  provided  in  geographic  coordinates,  and  are 
converted  to  the  sphere  trajectory  system  for  direct  comparison. 
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Table  4.1  Summary  of  Atmosphere  and  Wind  Profile  (AC13) 


AlTITUOt 

t-N  MIND 

M-6  MIND 

TIMP 

T8MP-U57# 

OtMSlTV 

DftHSnv  8 AT  10 

(KM) 

(M/S) 

(M/S) 

(OCC  K) 

(016  8) 

<*6/M**3) 

TO  US$476 

140 , 

-56.4 

34.1 

590.4 

30.7 

4 .6706-09 

1.193 

139. 

-58.3 

36.6 

684.5 

33.1 

4.666C-09 

1.164 

130. 

-67.9 

30.  1 

681  .3 

38.3 

6.1356-09 

1.168 

137. 

-60 .6 

3U.0 

677.9 

43.6 

8.4371-09 

1.183 

130. 

-54.4 

46.0 

674.1 

49.5 

6.7641-09 

1.140 

133. 

-53.3 

63.9 

669  8 

63.1 

6.12U-09 

1.120 

131. 

-53.3 

62.9 

665.2 

67.7 

6.6046-09 

1.103 

133. 

-52. U 

64.5 

6G0.8 

<2.8 

« . 9 1 16-09 

1.064 

133. 

-49.5 

57.3 

653. 6 

84.9 

7.3001 -09 

1  .070 

131. 

-44.9 

60.3 

545.6 

66.7 

7. 9186-09 

1.067 

130. 

-40.8 

60.4 

537.2 

•7.8 

9.6I0E-09 

1  .044 

139. 

-37.3 

57.4 

526.2 

67.0 

9.1976-09 

1.035 

130. 

-34.1 

52.8 

5 1  3 . 9 

64.9 

9.9b5t-09 

1  .020 

137. 

-30.8 

51 . 1 

503.  t 

64.5 

1.0836-00 

1.016 

130. 

-28.4 

65.0 

497.0 

69.0 

1.1656-06 

.995 

13o. 

-29.0 

61  .6 

496.1 

78.9 

1 .2416-00 

.962 

134. 

-30.7 

65.9 

499.3 

93. 1 

1.3126-08 

.019 

133. 

-35.6 

64.9 

502.3 

107.4 

1.3806-08 

.076 

133. 

-41.4 

61  .9 

502.5 

119.0 

■1 .4756-06 

.635 

131  . 

-47.2 

59.9 

497.2 

135.3 

1.587E-08 

.903 

130. 

-52.0 

67.1 

494.9 

124.9 

1 .7341-04 

.780 

1W. 

-57. 1 

68.7 

470.4 

•32.4 

1 .9096-08 

.761 

1  <0. 

-61 . 7 

67.  t 

450.6 

114.5 

3. 134L-D6 

.700 

1  17. 

-60. 0 

66.9 

428.8 

104.8 

2. 4106-08 

.742 

1 10. 

-89.9 

•3.3 

409.9 

97.9 

2.7196-04 

.fit 

1  10. 

-73.8 

43. U 

2»9 . 6 

89.6 

3.0986-08 

.722 

114. 

-73.3 

31.6 

lot  .8 

73.8 

3.631 6-08 

.730 

113. 

-70.  1 

19.0 

332.6 

66.5 

4.3326-00 

.746 

t  <3. 

•OJ.  3 

6.0 

3l»0.0 

44 . 9 

ft. 1616-09 

.766 

111. 

-49.0 

-14.4 

2u4 . 3 

42.  J 

9.0421-00 

.770 

110. 

-35.9 

-39.0 

362.9 

12.9 

7.6166-00 

.016 

too. 

-13.8 

-27.4 

335.1 

4.7 

4.7206-06 

.830 

10U. 

-15.0 

-8.4 

339.0 

15.7 

1 .096C-07 

.794 

107. 

3.0 

14.7 

207.9 

-9.9 

». 4551-07 

.00O 

too. 

3 »  .  to 

36.3 

178.3 

-34.8 

2.0146-07 

1.031 

10b. 

39.2 

24.0 

168.9 

-40.0 

2.8636-07 

1.103 

104. 

39.8 

10.9 

174.1 

-31.2 

3.0HC-07 

1.060 

103. 

35.9 

-5.2 

177.3 

-24.9 

3.5506-07 

1  .070 

103. 

50.3 

-16.6 

172.4 

-27.1 

4.4116-07 

1.121 

101. 

61.9 

-25.0 

174.0 

-23.2 

6.2906-07 

1.127 

too. 

61.0 

-28.7 

181.9 

-13.2 

6.0976-07 

1  .000 

99. 

62.9 

-39.3 

190.7 

-2.5 

0.9366-07 

1  .031 

96. 

49.4 

-32.1 

190.0 

-1.7 

8.2726-07 

t.626 

•7. 

63.9 

-39.6 

193.4 

-7.0 

1.0246-06 

1  .057 

90. 

72.9 

-43.6 

171.2 

-16. 1 

1 .3336-06 

1.I3M 

9b. 

97.6 

-29.8 

174.2 

-14.2 

t. 6766-00 

1 .131 

94. 

94.9 

-27.2 

196.8 

-1.0 

1.7696-06 

1.059 

93. 

78.1 

-23.3 

206.2 

19. 0 

1.6966-06 

.940 

93. 

72.2 

-29.3 

210.0 

23.0 

2.1666-06 

.913 

9f . 

63.4 

-29.0 

214.3 

27.6 

2.6048-06 

.070 

90. 

534 

-22.3 

21ft. i 

26.8 

2.8086*06 

.061 

99* 

60.1 

-34.6 

318.8 

32.7 

3. 3291-06 

.018 

Oh. 

63. 9 

-29.4 

221.1 

34.a 

8.8416-06 

.700 

97. 

44.2 

-35.7 

233.2 

46.3 

4.2161-06 

.734 

86. 

38.7 

-21.2 

226.6 

41.8 

4.9636-06 

.713 

es. 

38.0 

10.8 

226.4 

37.5 

6.6076-06 

.708 

•4. 

3.3 

35.3 

230.6 

39.8 

.6.5896-06 

.600 

93. 

13.9 

42.9 

727.® 

35.1 

7.7096-06 

.678 

83. 

18.6 

51.3 

2227 

26.0 

9.1436-06 

.601 

ttl  . 

11  .5 

51  .5 

226.8 

30.1 

1 .0426-06 

.661 

60. 

3.3 

49.0 

233.9 

36.3 

1 .1726-05 

.636 

79. 

6.  1 

97.0 

240.5 

40.Q 

1.3081 -06 

.606 

70. 

16.4 

ii.3 

342.2 

39.6 

1 . 4806-08 

.590 

77. 

39.9 

71.6 

240.7 

38.2 

1 .7206-09 

.804 

76. 

309 

43.1 

240.8 

34.1 

1 .6776-06 

.876 

7b. 

51.2 

46.6 

229.6 

30.1 

3.3971  05 

.601 

74. 

69.  1 

75.6 

220.7 

10.4 

2.8011-05 

.621 

73. 

64.5 

63.1 

331.3 
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3.3431-08 

.621 

73. 

63.3 

51  .6 

223.1 

7.8 

3.6726-05 

.621 

71. 

63.7 

40.6 

227.2 

10.4 

4, 3936*45 

.610 

70. 

72.5 

29.0 

230.2 

10.8 

6.0166-05 

.606 

69. 

98.6 

21.2 

232.2 

6.9 

6.7456-06 

.604 

60. 

69.6 

24.2 

231.7 

6.7 

8. 4*56-05 

.600 

67. 

58.  1 

13.3 

221  .6 

-6.0 

6. 0466-06 

.644 

6b. 

69.6 

3.4 

221  .3 

-6.3 

6.3086-05 

.668 

65. 

61 . 1 

-ft  .5 

225.  f 

-6.2 

1,07 1 (-04 

.667 

64. 

48.8 

-14.1 

227.7 

-8.4 

1.3286-04 

.660 

63. 

.3 

1.4 

249.0 

10.3 

1 .29*71-04 

.611 

63. 

2.0 

10.6 

254.3 

12.6 

1 .4466-04 

.601 

61 . 

9.2 

26.1 

254.0 

6.6 

1 ,6616-04 

.604 

60. 

7.4 

39.9 

252.8 

9.6 

1.8966-04 

.612 

6v. 

6.  7 

36.2 

263.6 

13.8 

2.0706-04 

.691 

66. 

7.9 

40.4 

272.3 

19. » 

2.2706-04 

.873 

67. 

8.8 

36.9 

274.1 

18.6 

2.8496-04 

.870 

96. 

6.8 

36.7 

271.8 

13.6 

2.8086-04 

.676 

96. 

9.8 

61.3 

183.4 

-7.4 

9.6388-94 

.633 
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In  this  manner  it  was  concluded  that  the  AC12  and  AC13  Z-accel- 
eration  values  were  about  15%  low,  and  the  revised  calibrations 
jointly  improved  in-track  wind  and  low  altitude  densities. 
Chemical  release  measurements  provide  good  wind  velocity 
estimates  at  the  higher  altitudes,  with  greatest  confidence 
placed  on  these  tests  under  conditions  of  uniform  winds, 
particularly  in  directions  transverse  to  the  lines  of  sight. 
Refinement  of  the  rocket  sphere  spin  rate,  suggesting  a  cross¬ 
track  axis  tilt,  would  be  possible  using  these  data. 


4.10  Evaluation  of  Reliability  of  Results 

Uncertainties  in  the  density  and  wind  solutions  arise  as 
discussed  above,  primarily  from 

a)  Calibration  coefficients 

b)  Trajectory  velocity  provided 
c>  Spin  axis  tilt  estimate 

d)  Spin  phase  orientation  estimate 


The  equations  that  have  been  described  in  Section  4 .6  for  the 
inversion  of  Eq.  4.4  are  also  directly  useable  for  determining 
the  sensitivity  of  the  solutions  to  variations  in  the  above  four 
parameters.  An  error  analysis  program  has  evolved  in  parallel 
with  the  complete  solution  algorithm,  and  turns  out  to  be  most 
useful,  not  only  to  establish  error-bars,  but  also  to  determine 
whether  any  combination  of  adjustments  of  the  parameters  would 
bring  results  closer  to  values  from  other  sources.  (Table  4.2) 
Further  advances  in  these  sensitivity  estimating  programs  are 
possible  in  order  to  include  other  variables  such  as  cross-track 
tilt,  and  for  implementation  of  alternate  solution  procedures  in 
the  80  km  transition  region. 


4.11  Auxiliary  Associated  Phenomena 


As  described  in  the  previous  sections,  studies  o-f  AC10,  AC12  and 

AC13  have  been  completed;  and  initial  processing  for  -flights  AC14 

thru  AC16  has  begun.  Results  and  procedures  as  now  implemented 

have  been  of  sufficient  resolution  quality  to  lead  to  the 

<2  3) 

reporting  of  basic  phenomena,  ’  the  identification  of  gravity 

<  4> 

waves,  and  the  association  of  density  gradient  reversals  with 
supersonic  auroral  fronts.^^  The  experiment  also  provides  a 
fertile  source  for  other  unscheduled  but  important  measurements 
in  the  lower  thermosphere,  such  as  estimating  the  drag  and  spin 
down  processes.  A  special  study  of  the  telemetry  AGC  signal  was 
conducted  for  the  sphere's  spinning  antenna,  in  order  to  fix  spin 
orientation.  It  was  shown  that  the  antenna  side-lobes 
introduced  unexpected  non-uniform  spin  modulation,  rendering 
phase  analysis  unreliable.  Future  daytime  flights  will  include 
a  sun-sensor,  and  accurate  spin  down  phase  determination  should 
be  possible. 

Atmospheric  turbulence  studies  were  also  performed  through 
calculation  of  Richardson  number  profiles  from  the  solved 
atmospheric  densi ty-temper ature  profiles.  This  parameter 
indicates  where  the  gradient  becomes  too  steep  to  sustain  a 
stable  density  profile. 
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5.0  Interactive  Graphics  Systems  for  AFGL  Applications 
5.1  Graphical  Functions 

Plotting  capabilities  very  similar  to  those  -for  conventional 
static  plotters  such  as  Calcomp  or  microfiche  are  valuable  to  a 
broad  range  of  investigators,  for  interactive  use  on  the  CDC-6600 
system  at  AFGL.  This  includes 

a)  Stacked  or  overlaid  X-Y  plotting  i.e.  more  than  one  plot 
di spl ayed . 

b)  Linear  or  logarithmic  axes  scales. 

c)  Symbol  and  connected  line  plotting  to  distinguish  the 
gr  aphs . 

d>  Additional  X-axes  for  concurrent  parameters. 

e)  Range  of  data  plotted  controlled  by  record  count  or 
dependency  variable  limits. 

f)  Flexible  heading  and  axis  labelling. 

g)  Running  average  and  other  simple  plot  smoothing  options. 

h)  Compatabi 1 i ty  between  off-line  and  on-line  data  base 
requirements  and  procedures. 

These  capabilities  have  been  implemented  for  conventional  X-Y 
plots  as  well  as  for  special  pictorial  presentations.  An 

(13) 

integrated  interactive  graphics  package  SUATEK  ’  incorporates 
the  above  and  other  useful  specialized  features.  Figures  5.1 
and  5.2  show  the  range  of  options  that  are  available  to  the 
SUATEK  user  in  menu  form.  Figure  5.1  lists  the  items  that  may 
be  set  by  the  user  in  order  to  select  variables  to  be  plotted, 
or  to  lay  out  the  plot  frame,  or  to  enable  options.  Activated 
options  may  work  in  extended  combinations,  and  Figure  5.2 
describes  each  one  briefly. 
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BRIEF  DESCRIPTION  OF  THE  TERNS 
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PIHN.PMX  -  HIN-NAX  ACCEPTANCE  RANGE  FOR  VARIABLE  IDP1 
BEFORE  CORRESPONDING  DATA  IS  PLOTTED 


DESCRIPTION  OF  OPTIONS  IPRRM . NCHK 
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Figure  5.2  Description  of  Activatable  SUATEK  Otions 


5.2  Analytical  Functions 


The  following  functions  can  be  performed  on  the  data  base  for 
direct  graphical  evaluation  and  revision: 

a)  Data  may  be  selectively  edited,  replaced,  and  filtered. 

b)  Data  may  be  plotted  conditionally,  depending  on  its 
range  or  that  of  some  other  significant  concurrent 
variable  (dependency  parameter) . 

c)  Data  may  be  sorted  according  to  one  or  more  variable 
prior  to  viewing. 

d)  Functional  combinations  of  variables  may  be  formed  as 
the  variable  to  be  plotted. 

e>  Polynomial  fits  to  segments  of  the  data  strearr  an  be 
applied  and  viewed. 

f)  Running  polynomial  fits  over  the  whole  data  s  _>  may 
be  applied  and  viewed. 

g)  Parameters  determined  during  any  of  the  above  analysis 
procedures  may  be  presented  for  recording  and 

eval uat i on . 

h)  New  data  formed  by  any  of  the  above  procedures  may 
replace  or  augment  the  original  data  base. 

5.3  Data  Base  Structure 

The  philosophy  for  a  basic  format  for  the  data  base  for  SUATEK 
has  been  the  capability  of  handling  data  streams  of  indeterminate 
length,  where  each  variable  to  be  considered  appears  in  each 
record,  in  standard  floating  binary  format.  The  only  other 
information  necessary  on  such  a  file  is  a  header  containing 

a)  Number  of  variables  per  record. 

b)  Alphanumeric  (axis  label)  per  variable. 

An  optional  plot  heading  label  may  be  added.  Figure  5.3  gives 
the  specification  for  SUATEK  files. 
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DESCRIPTION  OF  SUATEIC  FORMAT  FILES 
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This  specif i cat  ion  is  the  simplest  possible,  and  places  the 
•fewest  constraints  on  the  organization  and  development  of  new 
functional  options.  The  header  records  turn  out  to  be  self- 
descriptive  if  the  file  needs  to  be  examined.  An  important 
aspect  of  this  file  design  is  that  it  frequently  corresponds  to 
the  sequence  in  which  experimental  data  are  telemetered  and 
stored  on  tape.  A  variable  that  might  appear  say  once  in  64 
records  is  identified  by  another  subcom  frame  number  variable, 
which  is  necessarily  available  in  the  transmission.  Dependence 
plotting  with  the  appropriate  subcom  number  then  achieves  the 
desired  separation. 

The  on l y  additional  restrictions  placed  on  the  file  structure 
relate  to  two  special  values  for  each  variable.  Since  a 
variable  value  may  not  exist  in  concurrence  with  others  in  a 
record,  a  dummy  fill  or  IGNORE  value  is  required  to  which  the 
processing  must  be  completely  transparent.  In  other  cases  a 
variable  value  may  be  absent  or  consciously  removed,  and  where 
the  system  is  required  to  be  aware  of  this  break  in  the  data. 
Here  a  DELETE  value  is  required,  which  causes  a  break  while 
plotting,  and  which  may  be  automatically  replaced  during  certain 
curve  fitting  procedures.  The  numbers  assigned  to  DELETE  and 
IGNORE  are  -9.E99  and  -9.9E99  respectively,  values  which  are 
extremely  unlikely  to  be  encountered. 

The  main  objection  to  the  above  data  base  structure  may  be  the 
inefficient  use  of  disk  file  storage  when  a  few  bits  worth  of 
information  is  floated  and  stored  in  a  full  66-bit  CDC  word. 
SUATEK  could  be  modified  to  handle  either  the  conventional  or 
a  packed  file  format  automatically  through  a  flag  bit  in  the 
header.  A  packed  file  processing  package  available  in  the 
SUNY  library  could  be  adapted  for  this  implementation. 
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5.4  Implementation  of  the  Internal  File  Processing 


The  processing  functions  outlaid  in  Section  5.2  are  generally 
required  to  act  on  only  two  variables  <X,Y>,  even  though  other 
dependency  variables  may  determine  which  records  are  accepted  or 
rejected.  Subsequent  stages  for  any  complex  process,  e.g.  edit 
followed  bv  curve  fit,  are  handled  by  two  auxiliary  files  which 
contain  only  the  accepted  X,Y  values  and  their  record  number,  and 
which  are  used  solely  for  internal  reading,  processing,  and 
plotting  bv  SUATEK.  A  processing  pass  on  one  of  these  files 
creates  the  next  stage  data  on  the  other  file,  and  this  see-saw 
procedure  is  continued.  Details  of  each  analytical  procedure 
will  not  be  discussed,  as  they  are  extensions  of  conventional 
techniques.  During  a  pass  such  as  curve  fitting,  DELETE  records 
are  slipped,  additional  data  are  read  in  order  to  form  an 
adequate  window  for  the  fit,  and  the  DELETE  records  are  replaced 
to  form  a  continuous  plot.  Merging  of  the  final  reworked  data 
into  the  original  file  is  readily  accomplished  using  the  record 
number  saved  on  the  auxiliary  files. 


5.5  Implementation  of  the  Interactive  Modules 

acceptance  of  the  graphics  capabilities  by  researchers  for  their 
own  use  is  a  significant  consideration.  A  sel f -prompt i ng  plot 
specification  is  displayed  and  forms  the  initial  menu.  The  user 
may  enter  a  "?"  for  a  short  description  of  any  parameter  dis¬ 
played.  Plot  spec  modification  is  accomplished  by  entries  that 
are  nearly  identical  to  the  spec  display.  This  special  input 
section  has  the  virtue  of  permitting  gradual  build  up  of  each 
plot  specification,  and  of  readily  becoming  habitual  and 
unnoticed  by  the  typical  SUATEK  user.  Figure  5.4  is  a  summary 
guide  for  interactive  input.  This  scheme  has  also  been  adopted 
for  the  SUACON  contour  plotting  system.  (Section  5.6) 


r 
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SUATEK  INPUT  GUIDE 


tO  ENTER  OR  CHANGE  PARAHETER  VALUES  * 

HAKE  ENTRIES  EXACTLV  AS  DISPLAYED. 

FOR  EXAHPLEi  XHIN-8S, XHAX-160. *NYl-3  ETC. 

EHDEDDED  BLANKS  ARE  IGNORED 

SEPARATE  ENTRIES  BY  *,*  OR  BY  «CR** 

WITH  A  NAXINUH  OF  80  CHARS/LINE, 

AND  EACH  ENTRY  FULLY  CONTAINED  IN  A  LINE. 

TERH1NATE  ENTRY  BY  IN  ANY  POSITION* 

OR  BY  *1*  IN  POSITION  ONE. 


CONTROL  IS  AS  FOLLOWS t 

/  IN  POSITION  ONE  INITIATES  PLOTTING. 

/  IN  ANY  OTHER  POSITION  GIVES  A  SUHHARY 
•PLOT  SPEC*  STATUS  PRINT-OUT*  AND 
IN  POSITION  TWO  ALSO  PRESENTS  ALL  OPTIONS. 

1  IN  POSITION  ONE  INITIALIZES  PLOT  HIN-HAX  VALUES 
FRON  TAPE1  FOLLOWED  BY  *PLOT  SPEC*  STATUS  SUHNARV. 

U  ADDITIONALLY  DISPLAYS  ALL  VARIABLES  AND  OPTIONS. 

//  IN  POSITIONS  ONE-TWO  STOPS  THE  PROGRAN. 


AN  OPTION  IS  ACTIVATED  WHEN  SET  TO  NON-ZERO. 
(SUBSEQUENT  PROCEDURES  ARE  SELF-PRONPTING ) 


Figure  5«^  Summary  Guide  for  SUATEK  Interactive  Input 
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Output  prior  to  and  after  plotting  freely  uses  prompting  requests 
to  interact  with  the  user,  and  calls  for  specific  data  or  actions 
as  the  function  that  has  been  activated  progresses.  Cross  hairs 
adjustment  for  editing  proceeds  in  a  similar  manner,  with  the 
appearance  of  the  cross  hairs  being  the  signal  for  user  action. 

Realization  of  some  of  the  above  functions  requires  interchange¬ 
able  use  of  list  directed  and  formatted  print  statements,  as  well 
as  the  PLOTLIB  and  TEKSIM  libraries.  Compatibility  of  usage 
between  these  utility  packages  was  explored  with  the  systems 
group  at  AFGL,  and  a  calibration  pass  that  was  originally 
incorporated  is  now  not  required. 

Realization  of  any  of  the  interactive  graphics  packages  must  be 
through  the  CDC-6600  INTERCOM  system.  This  system  is  limited  to 
76 Kg  memory,  and  therefore  calls  for  an  Over l ay  or  Segmented 
organization  of  routines  that  need  not  be  co-resident.  This  is 
optimized  in  SUATEK  by  branching  into  one  of 
i)  the  input  section 

ii>  the  analytical  processing  sections  (which  may  have 
their  own  sub-branches) 
i 1 i >  the  plot  section 

Data  file  handling  and  common  blocks  are  resident  at  all  times. 
File  buffer  sizes  are  set  at  the  6 4  word  size  of  the  mass  storage 
PRU' s. 

The  NOS/BE  job  control  procedure  language  has  been  implemented, 
and  the  user  can  begin  to  plot  with  minimal  effort  —  viz.  one 
command:  BEGIN , SUA ,TEK ,M , HI SFI LE , I D=HI SNAME . 

where  M  is  the  Tektronix  unit  number  and  may  also  be  PEN  or  MIC 
for  off-line  plotting  purposes. 
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A*  mentioned  above,  pen-and-ink  or  microfiche  plots,  identical  to 
those  displayable  on  the  Tektronix,  may  be  initiated  using  SUATEK 
-files.  The  data  input  plot  requested  is  also  identical  to  the 
entries  that  are  keyed  in  on-line.  Thus  the  user  does  not  have 
to  learn  a  second  plotting  system.  Figure  5.5  owes  the  sample 
deck  set-up  -for  initiating  these  plots.  Four  cards  for  spec  1  a ' 
data  and  labelling  precede  the  normal  SUATEk  type  input.  Other 
di-fferences  are  that  the  off-line  <  SUAPEN)  version  does  not  allow 
interactive  input;  values  such  as  N  for  an  N-point  running 
average  must  be  input  where  the  prompt  would  normally  occur. 
Off-line  plotting  is  also  likely  to  consist  of  multiple  plots 
from  successive  files,  possibly  with  identical  plot  specifi¬ 
cations.  These  extra  facilities  are  implemented  in  SUAPEN. 

The  requirement  for  versatile  file  extraction  and  merging  capabi¬ 
lities,  arising  out  of  the  proliferation  of  SUATEK  format  files 
for  analyses  or  plots,  has  been  supported  specifically  by  program 
SUAMRG .  Specific  segments  and  variables  may  be  selected, 
relabelled,  and  merged  with  or  without  acceptance  of  special 
variable  control  values  such  as  DELETE.  Figure  5.6  summarizes 
SUAMRG  options,  and  shows  a  sample  execution. 


5.6  Three  Dimensional  Graphical  Systems 

Other  graphical  presentation  capabilities  are  valuable  for  multi¬ 
variable  files,  particularly 

i)  Contour  plotting 
ii)  Isometric  plotting 

Employing  different  approaches,  these  plot  forms  permit  introd¬ 
uction  of  another  dimension  for  viewing  and  interpreting  data. 
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JOBNM, CM7777? .  SUATEK  OFFLINE  PENPLOTS  ????  YOURNANE 
BEGIN, SUA.TEK, PEN, YOURSUFILE, ID«VOURID. 

EOF 

SCALE  FAC(DEF*1. )  1  IN  COL.  30  USES  PLOT  SPEC  SET  REPEATEDLY 

VOURNAHE  TEL  NO 

(BLANK  CARD) 

(BLANK  CARD) 

XL-180.,XHIN*240,XMAX-420  (PLOT  SPECS 

NV1-2,Y1MN-0.,V1MX-8.,YIL-1.  IN  STANDARD 

NV2«3,V2HN-0,Y2HX-.8,V2L-4,Y2OF-1.25  INTERACTIVE 

NV3-4,V3HN-0,V3HX-.8,V3L*4,Y3OF.S.S/  SUATEK  FORHAT) 


✓/ 

EOR 


(SIMILAR  ENTRIES  FOR 
FURTHER  PLOTS) 


Figure  5.5a  Sample  Deck  Set-up  for  Pen-and-ink  Plots 


PENPLOT  OR  MICROFICHE  PLOTS  HAY  BE  HADE  FROM 
SUATEK  FILES.  NO  NEU  INSTRUCTION  REPERTOIRE 
NEEDS  TO  BE  LEARNED,  SINCE  PLOT  SET-UP  INPUT 
IS  IDENTICAL  TO  THAT  FOR  ON-LINE  SUATEK. 


JOBNM, CH77777 .  SUATEK  MICROFICHE  PLOTS  7777  YOURNANE 

BEGIN, SUA,TEK, MIC, VOURSUFILE, ID* YOURID. 

EOR 

SCALE  FAC(DEF-1.)  1  IN  COL.  30  USES  PLOT  SPEC  SET  REPEATEDLY 

YOURNANE  TEL  NO 
MICROFICHE  TITLE1  (42  CHARACTERS) 

MICROFICHE  TITLES  (42  CHARACTERS) 

XL-10. ,XMIN-240,XMAX-260 
NY1*2,Y1HN*0.,Y1HX*8.,Y1L*1. 

NY2*3, Y2MN-0, V2MX* . 8 , V2L*  4 ,  Y20F*  1 . 25 
NY3*4,Y3HN*0,V3MX*.8,Y3L*4,Y3OF*S.5/ 

/ 

(SIMILAR  ENTRIES  FOR 
FURTHER  PLOTS) 


IN  COL.  50  FOR  MODULO  7  FRAMES 
(PLOT  SPECS 
IN  STANDARD 
INTERACTIVE 
SUATEK  FORHAT) 


// 

EOR 


Figure  5.5b  Sample  Deck  Set-up  for  Microfiche  Plots 
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SUAHRS  OPTIONS  FOR  HEROINS  TOPER  AND  TAPEB 

OPTION  1  —  INTERLEAVE  IASED  ON  MUTUAL  VARIABLE 
EQUIVALENCE  RECORDS  WITHIN  TOLERANCE 
OPTION  2  —  INTERLEAVE  BASED  ON  MUTUAL  VARIABLE 
INCLUDE  ALL  RECORDS 
OPTION  3  —  STACK  WRITE  TAPEA  RECORDS 
FOLLOWED  BY  TAPES  RECORDS 
OPTION  4  —  STACK  SEPARATED  BY  'DELETE'  RECORD 

SELECT  OPTION  (DEFAULT  0>1)  ***  0 


TAPEA  VARIABLES  FILTER  0  DO UN  LEG 

1  TINE 

2  ALTITUDE 

3  16 

4  FILTER 

5  COUNTS 

TAPEB  VARIABLES  FILTER  3  BOUN  LEG 
1  TINE 
ALTITUDE 
16 

FILTER 
COUNTS 

TO  MERGE  ENTER  NOVA,  NOVi  (0,0  TO  EXIT) 

ONE  COMMON  VARIABLE  (NOVA.NE.0  AND  NOVI.NE.f)  ALLOWED 

ENTER  NOVA,  NOVB  1^ 

ENTER  TOLERANCE  ANB  RANGE  (DEFAULTED, 0,i) 

ENTER  NOVA,  NOVI  2,0 
ENTER  NOVA,  NOVB  5.0 
ENTER  NOVA,  NOVB  0,5. 

ENTER  NOVA,  NOVB  0.0 

TOLERANCE  FOR  EQUIVALENCES  RECORDS  SET  AT  .15 

ENTER  ALPHANUMERIC  LABEL 
MERGED  R0D  AND  R3D  FILES 

REVISE  VARIABLE  LABELS?  <0*NO)  J_ 

REVISE  PLOT  MIN-MAX  RANGES?  <0«NO>  0. 

1  TINE 

2  ALTITUBE 

3  COUNTS 

4  COUNTS 

502  RECORIS  WRITTEN  TO  TAPE I 
END  SUAMRG 

Figure  5.6  SUAMRG  Options  and  Sample  Execution 


5.6.1  Contour  Plotting 

Characteristics  of  the  internet ationship  between  measurements  and 
causal  parameters  can  be  highlighted  through  these  presentations. 
NCAR  (Boulder,  CO)  has  developed  a  comprehensive  package  that 
offers  many  versions  of  these  capabilities  which  are  suitable  for 
off-line  use  but,  due  to  their  computer  memory  requirements  and 
large  range  of  options,  cannot  be  operated  in  on-line  (70Kg) 
mode.  Contour  plotting  is,  however,  a  presentation  that  offers 
valuable  insight  in  many  applications,  and  the  CONREC  routines  as 
adapted  by  Aiken  (AFGL) ,  were  implemented  for  on-line  Tektronix 
plotting.  This  system,  SUACON ,  allows  interactive  adjustment  of 
simple  options  such  as  frame  size,  axes  scales  and  labels,  and 
selection  of  contour  levels  to  be  plotted.  Figure  5.7(a),  (b) 
are  two  pages  from  the  current  user's  guide.  The  binary  data 
file  must  consist  of  a  header  and  successive  records  of  a  fixed 
number  of  X,Z  data  pairs: 


Record  1  -  MLBL(5) , L8L(3,2) 

where,  MLBL  is  a  50  character  alphanumeric  plot  label 
LBL  contains  30  char  X  and  Y  axes  l»b&ls 

Record  2  ...  NX, (X< I) ,Z< I) , 1=1 ,NX) 

where,  NX  =  number  of  X,Y  pairs  in  each  record 
X(I)  =  X  values,  must  be  equally  spaced 
2(1)  =  function  of  X,Y  being  contoured 
Successive  records  represent  equally  spaced  Y  values. 

A  maximum  of  50  records  with  NX  <  50  i s  presently  allowed. 

This  format  is  compatible  with  program  ISOTEK  which  produces 
isometric  plots  of  the  same  data. 
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8UACON  —  INTERACTIVE  CONTOURING 

Following  an  audit  o f  tho  input  data  file,  SUACON 
offer*  th#  user  a  number  of  default  or  activated!* 
option*  for  ponoratino  contour  plot*  of  th*  data. 

Th*  option*  and  th*  interaction  ar*  d**crib*d  using  a 
sample  session,  as  reproduced  on  the  following  oaoes. 


Audit  and  Initial  Presentation 

X  and  Z  min-max  value*  ar*  as  encountered  on  the  input  file. 

Y  value*  corresponding  to  each  record  do  not  occur  on  the  file 
and  may  therefore  be  redefined  during  th*  plotting  session. 
Default  Y  min-max  values  are  set  to  th*  X  min-max  values. 
Entering  C  directly  produces  the  default  plot  of  Example  1. 


Options  default  values  may  be  retained,  or  change*  may 
be  keyed  in  as  shown  underlined  in  th*  examples. 

XL,YL  specify  th*  plot  frame  i.e.  th*  length  of  th* 

X-  and  Y-  axes  in  inches  (default  ■  8  inches) 


XM1N ,XMAX  are  the  minimum  and  maximum  for  the  frame  range  and 
YMIN.YMAX  establish  the  plot  scale,  given  th*  axes  lengths 

XOFF.YOFF  ar*  th*  offsets  for  positioning  th*  plot  frame 


Options  below  are  off<»8>  but  may  be  selectively  activated(*l) 

ILABC  labels  the  contour  level* 

IDASHN  draws  negative  level  contours  dashed 

IMWN  labels  local  maxima-minima  levels 

ILBL  allows  modification  of  header  and  axes  labels 


IPICK.NCONT  ar*  used  for  Contour  Level  Specification  as  follows: 

Desired  contour  levels  to  be  drawn  must  always  be  identified, 

and  may  also  be  interactively  revised  after  viewing  the  plot. 

With  I  PICK  off,  NCONT  evenly  spaced  contour  levels  may  be 
automatically  selected  as  in  Example  Z. 

With  IPICK  on,  contour  levels  may  be  individually  specified 
as  in  Example  3, 

or .  previously'  set  contour  level*  may  be  listed 
and  selectively  revised  as  in  Example  4. 


Figure  5»7(a)  Excerpt  from  SUACON  User's  Guide 
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EXAMPLE  4.  ENTRIES  TO  PLOT  WITH  REVISED  CONTOUR  LEVELS 
(LEVELS  1  4  10  DELETED) 

***£ 

REPENT  PREVIOUS  PLOT?  (1-YES)J> 

REVISE  CONTOUR  LEVELS  NS  FOLLMIS 
ENCH  CURRENT  LEVEL  WILL  IE  DISPLAYED 
ENTER  -1  TO  DELETE  N  LEVEL 

•  TO  RETNXN  N  LEVEL 

1  TO  REVISE  THE  LEVEL 


LEVEL 

1  CONTOUR  VNLUE 

•9S0SE+M 

7-1 

LEVEL 

2  CONTOUR  VNLUE 

.steoE+ee 

T 

LEVEL 

3  CONTOUR  VNLUE 

.7S0SE+M 

4 

LEVEL 

4  CONTOUR  VALUE 

.SIME4H 

?• 

LEVEL 

S  CONTOUR  VALUE 

.2SHE4M 

?I 

LEVEL 

6  CONTOUR  VALUE 

%25NEtM 

?I 

LEVEL 

?  CONTOUR  VNLUE 

-.SOOOE+OO 

?± 

LEVEL 

g  CONTOUR  VNLUE 

-.7SS0E+0S 

?A 

LEVEL 

9  CONTOUR  UNLUE 

?£ 

LEVEL 

IS  CONTOUR  UNLUE 

-.KIKtH 

7-1 

DO  YOU 

UISH  TO  ADD  NNV  LEVELS?  C  1-YES )?JL 

contour  nor  test 

I.Of  -  »-  —4 - 4- - 4--  >  «  ■■  «  I  > 


i.u 

X-COORO 


Figure  5»7(b)  Excerpt  from  SUACON  User's  Guide 


The  Procedure 


After  ETL,???  for  extended  interactive  sessions,  SUaCON  may  be 
set  up  and  activated  for  Tektronix  and  associated  Pen  and  Ink  or 
Microfiche  plotting  with  a  single  command  as  follows: 

BEGIN , CON , TOUR ,MODE , pf  n , 1 D=pf i d . 

where  MODE  =  1  or  2  (Tektronix  14)  for  Interactive  use, 

PEN  or  MIC  (off-line  not  implemented) 

If  ID  is  omitted  a  local  file  "pfn"  is  assumed, 
and  a  local  file  TAPE1  is  created. 

If  pfn  is  also  omitted,  the  default  lfn  is  TAPE1. 


5.6.2  Isometric  Plotting 

Another  package  (ISOTEK)  developed  from  basic  principles  of 

hidden  line  removal,  and  using  an  identical  data  base  format  to 

the  one  for  SUACON,  was  developed  for  isometric  presentations  of 

three-dimensional  data.  The  system  simplifies  the  basic 
( 2) 

approach  of  Wray  by  allowing  interactive  selection  of  block 
min-max  ranges  and  therefore  the  viewing  orientation  and  scale. 
Block  depiction,  tick  marks,  cross-hatching,  and  other  simple 
options  are  implemented.  Figure  5.8  shows  sample  usage. 
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commas-  iaa 

ENTER  STMT  ANO  END  RECORD  NUMBERS  (0EF-«.*>  >*,» 

ENTER  RIOT  I  LOCK  UNITS  (BEf-R.R) 

XAIN.XNAX  >-2.S,7.S 
YN1N.YNAX  >-4.1 8 

TO  HIDE  LINES  INTERPOLATE  UlTH  EUENLV  SPACES  POINTS 
ENTER  NUNSER  Of  POINTS  PER  LINE  (O-PLOT  RAW  SATA) 

ENTER  NUNSER  Of  CROSS  LINES  WANTED  (DCF-*)  >51 

S*  GRID  LINES  HADE  SEEN  COPIED 

UlTH  SUATEK  FORMAT  ON  TAPE  I 
STOP 

******  MAXIMUM  EXECUTION  FI. 

7.7*3  CP  SECONDS  EXECUTION  TINE. 

COMMAS-  SUATEK 2 


>25* 


x 


COMMAS-  ISOTEK 

ENTER  START  AND  END  RECORD  NUNiERS  (D£F-«.«) 


ENTER  PLOT  BLOCK  LIMITS  <D£F-*.*> 
XHIN.XMAX  >-S-S.P.S 
VNIN. VNAX  >-1.6 


>2S.4B 


TO  HIDE  LINES  INTERPOLATE  WITH  EUENIV  SPACES  POINTS 
ENTER  NUMBER  OF  POINTS  PER  LINE  <*-PLOT  RAW  DATA) 

ENTER  NUNSER  OF  CROSS  LINES  WANTED  (DCF-*)  >£ 

24  GRID  LINES  HAVE  SEEN  COPIED 
UlTH  SUATEK  FORMAT  ON  TAPE  1 
STOP 

060*0*  MAXIMUM  EXECUTION  FL. 

3.368  CP  SECONDS  EXECUTION  TINE. 

COMMAND-  SUATEK2 


>2** 


Figure  5*8  ISOTEK  -  Interactive  Isometric  Sample  Plots 
(An  intermediate  file  is  created  for  SUATEK) 
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